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ABSTRACT 

Context. NGC 2264 is a young stellar cluster (~ 3 Myr) with hundreds of low-mass accreting stars that allow a detailed analysis of 
the accretion process taking place in the pre-main sequence. 

Aims. Our goal is to relate the photometric and spectroscopic variability of classical T Tauri stars to the physical processes acting in 
the stellar and circumstellar environment, within a few stellar radii from the star. 

Methods. NGC 2264 was the target of a multiwavelength observational campaign with CoRoT, MOST, Spitzer, and Chandra satellites 
and photometric and spectroscopic observations from the ground. We classified the CoRoT light curves of accreting systems according 
to their morphology and compared our classification to several accretion diagnostics and disk parameters. 

Results. The morphology of the CoRoT light curve reflects the evolution of the accretion process and of the inner disk region. 
Accretion burst stars present high mass-accretion rates and optically thick inner disks. AA Tau-like systems, whose light curves are 
dominated by circumstellar dust obscuration, show intermediate mass-accretion rates and are located in the transition of thick to 
anemic disks. Classical T Tauri stars with spot-like light curves correspond mostly to systems with a low mass-accretion rate and 
low mid-IR excess. About 30% of the classical T Tauri stars observed in the 2008 and 2011 CoRoT runs changed their light-curve 
morphology. Transitions from AA Tau-like and spot-like to aperiodic light curves and vice versa were common. The analysis of the 
Her emission line variability of 58 accreting stars showed that 8 presented a periodicity that in a few cases was coincident with the 
photometric period. The blue and red wings of the Her line profiles often do not correlate with each other, indicating that they are 
strongly influenced by different physical processes. Classical T Tauri stars have a dynamic stellar and circumstellar environment that 
can be explained by magneto spheric accretion and outflow models, including variations from stable to unstable accretion regimes on 
timescales of a few years. 
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1. Introduction 

Classical T Tauri stars (CTTSs) are young, low-mass stars (M=, < 
2 Mo), with spectral types from F to M. They are surrounded by 
a circumstellar disk from which they are still accreting mate¬ 
rial. They present strong and broad emission lines in their spec¬ 
tra and show emission excess with respect to the stellar photo¬ 
sphere that goes from the radio to the ultraviolet (Bouvier et al. 
2007b). CTTSs have strong magnetic fields (Johns-Krull et al. 
1999; Johnstone et al. 2014) that disrupt the accretion disk at a 
few stellar radii from the star and channel the accreting material, 


^ Tables 2 and 3 are only available in electronic form at the 
CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via 
http://cdsweb.u-strasbg.fr/cgi-bin/qcat7J/A-hA/ 


forming accretion columns. The accreting gas hits the stellar sur¬ 
face and creates hot spots. CTTSs also present cold spots at the 
stellar surface; these are caused by magnetic activity (Bouvier 
et al. 1995). Part of the gas in the inner disk region is ejected 
as a wind from the star-disk system along open magnetic field 
lines that may form collimated jets (e.g., Ferreira et al. 2006). In 
a few million years, before reaching the main sequence, CTTSs 
lose their disks and become weak-lined T Tauri stars (WTTSs), 
which no longer show detectable signs of accretion (e.g., Meyer 
2009). 

A characteristic of CTTSs is the photometric and spectro¬ 
scopic variability at various wavelengths. The photometric vari¬ 
ability occurs from X-ray to infrared on a timescale from a few 
minutes to several years and is usually irregular (e.g., Appen- 
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zeller & Mundt 1989). Some stars, however, show periodic be¬ 
havior, which may be due to the presence of stable cold and 
hot spots on the stellar surface or to circumstellar dust extinc¬ 
tion, as observed in the classical T Tauri star AA Tau (Bouvier 
et al. 2007a; Alencar et al. 2010). The analysis of the photo¬ 
metric variability of CTTSs allows the determination of cold or 
hot spot characteristics and the typical timescale of the physical 
processes that cause each type of variability (dynamo, accretion, 
star-disk interaction). We can also estimate the line-of-sight dust 
distribution in the inner disk in favorable star-disk inclinations, 
when the inner disk occults the star as the system rotates. Spec¬ 
troscopic variations are also present, and emission lines can vary 
in shape and intensity on a timescale of hours to days (e.g., Johns 
& Basri 1995a; Muzerolle et al. 1998; Costigan et al. 2014). The 
study of the spectroscopic variability of CTTSs can be related to 
the predictions of magnetospheric accretion models and magne- 
tohydrodynamical simulations, which may include the accretion 
and the wind components, as well as the star-disk interaction. 

CTTS also show infrared excess emission that indicates the 
presence of a circumstellar disk (e.g., Teixeira et al. 2012). This 
emission can be used to estimate the amount of dust in the sys¬ 
tem and relate disk and accretion evolution. The more evolved 
the star-disk system, the smaller the amount of circumstellar ma¬ 
terial available, and consequently, the lower the emission excess 
in the infrared. 

Magnetospheric accretion is the standard model to describe 
the accretion process in CTTS (Shu et al. 1994; Hartmann et al. 
1994). Magnetohydrodynamic simulations predict that magne¬ 
tospheric accretion can occur in stable and unstable regimes 
(Kurosawa et al. 2008; Kulkarni & Romanova 2008,2009; Kuro¬ 
sawa & Romanova 2013). In the stable regime, accretion occurs 
through two main accretion funnels (one in each hemisphere), 
and periodic spectroscopic and photometric variations are ex¬ 
pected (Kurosawa & Romanova 2013), since there is a global 
organization of the accretion geometry. In the unstable regime, 
several accretion streams are formed that appear at random loca¬ 
tions, creating multiple hot spots on the surface of the star. This 
accretion regime can be maintained by the Rayleigh-Taylor in¬ 
stability, which acts at the magnetosphere-disk interface (Kuro¬ 
sawa & Romanova 2013) and causes stochastic photometric and 
spectroscopic variability. If the Rayleigh-Taylor instability is 
weak, both the unstable and stable accretion regimes may coex¬ 
ist. Accretion will then occur mainly through two accretion fun¬ 
nels, but random accretion funnels can also appear in the mag¬ 
netosphere (Kulkarni & Romanova 2008). These predictions can 
be checked with photometric and spectroscopic observations of 
CTTSs that span different timescales, from days to years, as the 
ones discussed in this work. 

We analyze young stars belonging to NGC 2264, a young 
stellar cluster (< 3Myr) located at a distance of ~ 760 pc from 
the Sun (Sung et al. 1997). This cluster shows evidence of on¬ 
going star formation, such as the presence of molecular out¬ 
flows and Herbig-Haro objects. In the pioneering work of Herbig 
(1954), 84 pre-main sequence stars were found in NGC 2264 
with Ha emission and were classified as T Tauri stars. Since 
then, and because of its proximity and low extinction in our line 
of sight, NGC 2264 has been the subject of many observational 
campaigns, from the radio to X-rays (e.g.. Rebull et al. 2002; 
Lamm et al. 2004; Dahm & Simon 2005; Dahm 2008; Teixeira 
et al. 2012). About 1000 pre-main sequence stars have already 
been confirmed as members of NGC 2264 (Sung et al. 2009). 

NGC 2264 was observed with the Convection Rotation and 
planetary Transits (CoRoT) satellite during 23 consecutive days 
in 2008 (from March 7 to 30), and photometric data with high 


cadence and high signal-to-noise ratio were obtained for ~ 300 
known cluster members. The analysis of these data yielded stud¬ 
ies in different research areas such as astroseismology (Zwintz 
et al. 2013), stellar rotation (Affer et al. 2013), binary systems 
(Gillen et al. 2014), and accretion (Alencar et al. 2010). In the 
latter, it was shown that the light curve exhibited by the clas¬ 
sical T Tauri star AA Tau, which is due to obscuration of the 
photosphere by circumstellar material in the inner disk region, is 
common among the cluster members. This opens the possibility 
of studying the inner disk evolution with photometric variabil¬ 
ity at different wavelengths. It was also shown that the CoRoT 
light curve morphology is related to the inner disk evolution and 
consequently to the accretion process. 

A second observational campaign of NGC 2264 was orga¬ 
nized in 2011 (Cody et al. 2013), including optical, infrared (IR), 
and X-ray simultaneous observations of the accreting and non¬ 
accreting members of the cluster. We present this in the next 
section. Using this new data set, we discuss the observed pho¬ 
tometric and spectroscopic variability and its relation to magne¬ 
tospheric accretion model predictions. We classify the CoRoT 
light curves of accreting systems according to their morphology 
and compare this classification with accretion diagnostics such 
as Ha and Hel 6678A emission from FLAMES at the Very Large 
Telescope (VLT), UV excess from MegaCam at the Canada- 
France-Hawaii Telescope (CFHT), and disk parameters such as 
IR excess. A set of non-accreting members of NGC 2264 is used 
as a control sample to define the role of accretion in the observed 
correlations. Results from previous observations of NGC 2264 
with the CoRoT satellite in 2008 are compared with the 2011 
campaign to analyze the dynamical nature of accretion and the 
inner disk evolution. We analyze the variability in light-curve 
morphology between the two CoRoT campaigns and relate it to 
the proposed stable and unstable accretion scenarios in the lit¬ 
erature. Cody et al. (2014) presented a detailed description of 
the photometric data sets and morphological classification of the 
various types of CTTS light curves observed with CoRoT and 
Spitzer in the 2011 campaign. 

The paper is organized as follows. In Sect. 2, we present the 
2011 observational data of NGC 2264 and the reduction pro¬ 
cedures. The CTTS sample selection criteria are presented in 
Sect. 3. In Sect. 4 we discuss the morphological classification of 
the CoRoT light curves, following the classification scheme pro¬ 
posed by Alencar et al. (2010). In Sect. 5 we compare the CoRoT 
light-curve classification to the CFHT photometric analysis un¬ 
dertaken by Venuti et al. (2014). In Sect. 6 we analyze the accre¬ 
tion rate obtained directly from the Ha line flux. In Sect. 7 we 
morphologically classify the Ha line, as proposed by Reipurth 
et al. (1996). In Sects. 8 and 9 we obtain the periodicities of the 
Ha and Hel 6678 A lines. We discuss the presence of unstable 
and stable accretion regimes in Sect. 10, and we analyze the cor¬ 
relation between the different Ha line structures (emissions and 
absorptions) in Sect. 11. In Sect. 12 we discuss some interesting 
objects to be analyzed individually in a future work. Section 13 
presents our conclusions. 


2. Data and reduction 

NGC 2264 was observed with the CoRoT satellite during 40 
days, from December 2, 2011 to January 10, 2012, providing 
exquisite photometry for ~ 500 probable cluster members. Si¬ 
multaneously with the 2011 CoRoT run, the cluster was ob¬ 
served with the Spitzer satellite at 3.6 and 4.5 pm for 30 days and 
for 3.5 days with the Chandra satellite. This campaign was called 


Article number, page 2 of 33page.33 



A. R Sousa et aL: CSI2264: Accretion process in classical T Tauri stars in the young cluster NGC 2264 


Coordinated Synoptic Investigation of NGC 2264 (CSI 2264)^ , 
and details about the data acquisition and reduction can be found 
in Cody et al. (2014). We also obtained u- and r -band observa¬ 
tions of the cluster with MegaCam/CFHT from February 14 to 
28, 2012. They were combined with observations of a first Mega- 
Cam campaign in 2010 to yield photometry in the ugri bands 
for a large sample of stars in the NGC 2264 region. A detailed 
description of the CFHT observations is given in Venuti et al. 
(2014). 

As a complement to the CSI 2264 campaign, we obtained 
FLAMES spectra at VLT/ESO of 58 CTTSs and 34 WTTSs. 
These stars were selected based on the CoRoT light-curve clas¬ 
sification of 2008 (Alencar et al. 2010), to include most of 
the A A Tau and spot-like systems. FLAMES is a high- and 
medium-resolution multi-object spectrograph, covering a field 
of view of 25’ in diameter. We obtained 20 to 22 spectra of 
each target in two different fields, distributed in 20 days from 
December 2011 to February 2012, part of them simultaneously 
with the CoRoT observations. We used the HR15N setup for 
the FLAMES/GIRAFFE spectrograph, centered at 6650 A and 
covering the 6470 A < A < 6790 A spectral region with a 
resolution of R = 17,000. This spectral region includes Ha, 
and the HeI6678A spectral lines that are analyzed in detail 
in the next sections. The observed spectra were reduced using 
the GAS GANG reduction package from ESO. We performed 
the standard procedure of bias subtraction, fiat-field correction, 
wavelength calibration, and barycentric velocity correction for 
all spectra and later normalized their continuum. In Table 1, we 
summarize the observations of the CSI 2264 campaign that are 
discussed in the following sections. 

NGC 2264 is located in front of an emission nebula that 
is seen spectroscopically as a series of emission lines, such 
as [Nil], [SII], and Ha, superposed on the stellar spectra. The 
Ha nebular contribution is very narrow and in general much 
smaller than the accretion contribution to the Ha emission of 
most CTTSs, presenting a mean equivalent width of the Ha neb¬ 
ular contribution of 10.1 + 0.5 A, averaged over all the sky fibers. 
However, in a few spectra, the nebular emission is completely 
blended with the emission from the star, making it impossible 
to identify whether a component is of nebular or accretion ori¬ 
gin. This sometimes makes removing the nebular emission from 
the total spectrum challenging. In general, this is a difficult task 
when using fiber spectrographs, as is the case of FLAMES. The 
nebular contribution is quite variable across the sky, and the sky 
fiber closest to the stellar fiber does not always present the same 
nebular contribution as the stellar fiber. Another difficulty is that 
some targets are very weak and present a spectrum comparable 
in intensity to the sky contribution. 

To remove the nebular contribution from the stellar spectra, 
we subtracted the sky spectrum of the nearest sky fiber from the 
stellar spectrum of each star (see Fig. 1). The subtraction of the 
nebular emission is sometimes faulty, and a residual nebular line 
often remains. Even when the sky subtraction leaves residuals 
in the stellar spectrum, the contamination from the nebula is 
normally distinguishable from the stellar emission because the 
sky emission is narrower and weaker than the emission of most 
CTTSs. 


^ All of the CoRoT and Spitzer light curves can be viewed and down¬ 
loaded from http://irsa.ipac.caltech.edu/data/SPITZER/CSI2264/ 



Fig. 1. Her emission lines of a CTTS before (left) and after (right) sky 
subtraction that show the removal of the nebular emission contribution. 
The nebular component is the very sharp, narrow emission seen near 
the line center in the left panel. Different colors correspond to different 
observing nights. 


3. CTTS and WTTS samples 

We selected the CTTSs in NGC 2264 from the known cluster 
members according to photometric and spectroscopic criteria. 
We measured the Ha equivalent width (EWh^) and Ha width 
at 10% of maximum intensity (W10%HQr) in the FLAMES spec¬ 
tra after removing the nebular contribution. We used published 
values of these parameters for the stars we did not observe spec¬ 
troscopically, as shown in Table 2. To test the influence of the 
nebular emission on the measured values of EWh^, we calcu¬ 
lated the Ha equivalent width before and after subtracting the 
nebular part, from which we obtained a mean difference of about 
10%. This difference does not affect our CTTS classification. 

According to White & Basri (2003), CTTS have an EWhot 
higher than a threshold that depends on the stellar spectral type. 
A young star is considered a CTTS if it presents an EWhqt higher 
than 3 A for K0-K5, 10 A for K7-M2.5, 20 A for M3-M5.5, 
and 40 A for M6-M7.5. Stars that show W10%Ha higher than 
270 km/s are also considered CTTSs, following White & Basri 
(2003). Fang et al. (2009) proposed more stringent limits to sep¬ 
arate WTTSs from CTTSs that extend to late-M spectral type. 
Since we do not have stars in this spectral type region, we de¬ 
cided to use the criteria proposed by White & Basri (2003) to 
select CTTS. We selected 220 CTTSs in NGC 2264 that satisfy 
at least one of the above spectroscopic criteria. 

Using UV excess obtained with MegaCam/CFHT, Venuti 
et al. (2014) found 66 new CTTS candidates in NGC 2264. 
These stars do not have FLAMES spectra and were not previ¬ 
ously classified as CTTS in the literature or are unknown NGC 
2264 members, but since they present considerable UV excess 
with respect to the WTTS locus, they were added to our CTTS 
list. This means that based on photometric and spectroscopic ac¬ 
cretion criteria, 286 stars are thought to be CTTSs in NGC 2264. 

Mon-000056 and Mon-000358 were initially not included 
in our CTTS list, but they show CoRoT light curves in 2008 
and/or 2011 that resemble the AA Tau light curve (see Sect. 
4). We would therefore expect them to be CTTSs, but they 
were not in the fields of our FLAMES observations and show 
little UV excess. These stars were nevertheless observed with 
FLAMES/GIRAFFE as part of the Gaia/ESO public spectro¬ 
scopic survey (Randich et al. 2013; Gilmore et al. 2012). We 
downloaded the spectra from the ESO science archive and re¬ 
duced them following the same steps as listed in Sect. 2. The 
Gaia/ESO spectra show that Mon-000056 presents EWh^ = 
2.9 A and W10%HQr = 276 km/s, while Mon-000358 shows 
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Table 1. Observations from the Coordinated Synoptic Investigation of NGC 2264 


Telescope 

Instrument 

Dates 

Bands 

CoRoT 

E2 CCD 

2011-Dec-1 to 2012-Jan-3 

3000- 10000 A 

VLT 

FLAMES/GIRAFFE 

2011-Dec-4 to 2012-Feb-29 

6470-6790 A 

CFHT 

MegaCam 

2012-Feb-14 to 2012-Feb-28 

u and r 



jKcrrs WTTS 


100 r 


^ 10 r 
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Fig. 2. Spectroscopic criteria used to select CTTSs and WTTSs in NGC 
2264. CTTSs are shown as asterisks and WTTSs as plus signs. Top: Her 
equivalent width vs. spectral type for our sample of stars. Dashed hori¬ 
zontal lines represent the White & Basri (2003) classification threshold 
and solid horizontal lines the Fang et al. (2009) proposed criteria. Bot¬ 
tom: Her equivalent width vs. Her width at 10% of maximum intensity. 
The dashed vertical line is drawn at 270 km/s, the separation between 
CTTSs and WTTSs, according to White & Basri (2003). A few stars do 
not satisfy one spectroscopic CTTS criterion, but they were classified as 
CTTSs either based on the other spectroscopic criterion or on the UV 
excess. 


EWhot = 32.8 A and W10%Ha = 395 km/s, which allows us 
to classify them as CTTSs. Our final CTTS list therefore con¬ 
tains 288 stars; it is shown in Table 2. Of these, 58 stars were 
observed by FLAMES and 157 were observed by CoRoT in the 
2008 and/or 2011 campaign, including 12 CTTSs observed by 
CoRoT only in 2008, 63 only in 2011, and 84 CTTSs that were 
observed in both campaigns. 

We also selected 431 WTTSs among the confirmed members 
of NGC 2264 that we used as a comparison sample to the accret¬ 
ing CTTS. We added 257 WTTSs to the list based on the Her 
equivalent width or W10%HQr criteria of White & Basri (2003), 
previously discussed, and 174 WTTSs were selected based on 
the criteria discussed in Venuti et al. (2014). Information about 
the WTTS sample can be found in Table 3. Of the 431 WTTSs, 
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34 were observed by FLAMES and 308 stars were observed by 
CoRoT in the 2008 and/or 2011 campaign. We did not analyze 
the CoRoT data for WTTSs; that is beyond the scope of this pa¬ 
per. 

In Fig. 2, we plot the spectroscopic criteria used to dis¬ 
tinguish CTTSs from WTTSs, as proposed by White & Basri 
(2003) and Fang et al. (2009), together with the Ha data for the 
58 CTTS and 34 WTTS observed with FLAMES. The spectro¬ 
scopic data of the stars that were not observed with FLAMES 
come from Dahm & Simon (2005). A few stars do not satisfy the 
EWhot CTTS criteria, but they were classified as CTTSs based on 
W10%HQr (Eig. 2, bottom panel) and/or UV excess. 

4. Morphology of the CoRoT light curves 

The 2008 CTTS light curves were previously classified morpho¬ 
logically as spot-like, AA Tau-like, or non-periodic by Alen- 
car et al. (2010). The spot-like light curves display sinusoidal 
variations that are due to large, cold spots at the stellar surface, 
which are very stable on the timescale of the observations (about 
three weeks). The AA Tau-like light curves show a well -defined 
maximum interrupted by periodic minima that vary in width and 
depth from cycle to cycle, like AA Tau itself. They correspond 
to systems whose light-curve variability is caused by a periodic 
occultation of the stellar photosphere by circumstehar dust (see 
McGinnis et al. 2015; Alencar et al. 2010, for more informa¬ 
tion about A A Tau-like stars in NGC 2264 and Eonseca et al. 
(2015) for a specific AA Tau-like system, V354 Mon - Mon- 
000660). Non-periodic light curves may be due, for example, to 
accretion bursts (Stauffer et al. 2014) or random circumstehar 
dust obscuration. Some non-periodic light curves are difficult to 
assign to a major physical phenomenon, however, and are prob¬ 
ably the result of different superposed variability processes, such 
as variable accretion, dust obscuration, and spot variations. We 
show below by comparing our light-curve classification with ac¬ 
cretion diagnostics and inner disk parameters that the morpho¬ 
logical classification of the CoRoT light curves is related to the 
evolution of the accretion process and of the inner disk region. 

We morphologically classified the 145 CTTSs light curves 
observed in 2011 by CoRoT, and in Fig. 3 we show examples 
of these light curves. In Fig. 4, we present the periodic light 
curves of Fig. 3 folded in phase with the period obtained with the 
Scargle periodogram, as modified by Home & Baliunas (1986). 
The stability of the spot-like light curve modulation is evident, 
which contrasts with the high variability of the typical AA Tau- 
like light curves. A detailed statistical and morphological anal¬ 
ysis of the light curves of stars that present infrared excess ob¬ 
served in our campaign with CoRoT and Spitzer are discussed 
in Cody et al. (2014), where the light curves are separated into 
twelve different classes that represent diverse physical processes 
and geometric effects. We list in the last column of Table 2 the 
Q parameter defined by Cody et al. (2014) that can be used as a 
metric to distinguish periodic and quasi-periodic (2 < 0.6) from 
aperiodic (2 > 0.6) light curves. 

Table 4 shows some morphological statistics of the CTTS 
light curves observed with CoRoT in 2008 and 2011. The num- 
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Mon-000103 Mon-000017 



JD-2450000 JD-2450000 


(a) 


(b) 


Mon-001140 



5890 5900 5910 5920 5930 5940 
JD-2450000 


Mon-001296 



5890 5900 5910 5920 5930 5940 
JD-2450000 


(C) 


(d) 


Mon-000676 



5890 5900 5910 5920 5930 5940 
JD-2450000 


Mon-000567 



5890 5900 5910 5920 5930 5940 
JD-2450000 


(e) 


(f) 


Fig. 3. CoRoT light curves of CTTSs from 2011 showing different 
photometric variabilities. Panels (a) and (b) correspond to spot-like light 
curves, (c) and (d) to AA Tau-like light curves, (e) is a non-periodic 
light curve caused by obscuration by circumstellar material, and (f) is a 
non-periodic light curve dominated by accretion bursts. The light-curve 
magnitude was calibrated using the R filter with a zero point of 26.74 
mag, as described in Cody et al. (2014). 


ber of periodic CTTSs did not vary considerably from one run to 
the other. However, the number of CTTSs observed with CoRoT 
increased in 2011, and more non-periodic stars were observed in 
2011 than in 2008. This was expected, since most of the known 
periodic variable systems were already included in the 2008 run. 
The number of AA Tau-like systems identified in 2008 differs 
from that of Alencar et al. (2010), since we were able here to 
use the simultaneous FLAMES spectroscopy to refine our clas¬ 
sification. In addition, Stauffer et al. (2015) identified a new sub¬ 
class of stars with variable extinction light curves that show peri¬ 
odic, shallow, and short-duration flux dips that are approximately 
Gaussian in shape. There is some overlap between this new class 
and the classical AA Tau systems, since a few AA Tau-like stars 
also present narrow dips superposed on the broad deep minima, 




Fig. 4. Periodic CoRoT light curves from Fig. 3 folded in phase. The 
spot-like light curves present stable variations (top) that contrast with 
the high variability of the AA Tau-like light curves (bottom). Different 
colors represent different rotational cycles. 


as explained in their paper. Light curves labeled “unclassified” 
in Table 4 correspond to the light curves that we were unable to 
clearly associate with one major physical phenomenon. 

Table 4. Morphological classification of CTTS light curves observed 
with CoRoT in 2008 and 2011. 


Light curve type 

N® of stars-2008 

N® of stars-2011 

Periodic or quasi-periodic 

Spot 

24 

19 

AA Tau 

15 

18 

Eclipsing binary 

2 

2 

Unclassified 

0 

4 

Non-periodic 

Random extinction 

13 

13 

Accretion bursts 

12 

20 

Unclassified 

30 

69 

Total 

96 

145 


4.1. Morphological changes from 2008 to 2011 

Several stars exhibited the same light-curve morphology in 2008 
and 2011, but ~ 30% of the stars presented variations in the 
morphology of their light curves. We show in Table 5 the varia¬ 
tion of morphological classification of the 84 CTTS light curves 
that were observed in both campaigns. The transition from AA 
Tau to non-periodic classification and vice versa is common, but 
we did not observe a transition from AA Tau to spot-like or 
from spot-like to A A Tau. A few spot-like systems became ir¬ 
regular, and the opposite also occurred. In Fig. 5 (a) and (b), 
we show the 2008 and 2011 light curves of the CTTS Mon- 
000928, whose light-curve classification changed between the 
two runs. In 2008, Mon-000928 was classified as an A A Tau-like 
system with a maximum brightness periodically interrupted by 
minima, while in 2011 it was clearly irregular with no detected 
periodicity. These variations show that the inner disk structure 
that is responsible for the stellar occultation can move from a 
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well-organized geometry with a stable inner disk warp to a ran¬ 
dom dust distribution in just a few years. These results may also 
hint at the typical timescale of stability of the accretion regime, 
where the accretion geometry evolves from a main accretion fun¬ 
nel in each hemisphere, the base of which may correspond to the 
stable inner disk warp to an unstable accretion scenario, where 
random accretion funnels are formed (Kurosawa & Romanova 
2013; McGinnis et al. 2015). 

In young low-mass stars, cold spots often occupy a substan¬ 
tial fraction of the stellar surface and may last for months to years 
(e.g., Herbst 1989; Vogt et al. 1999). If the spots indeed did not 
evolve significantly over a period of few years, then the light 
curves would also be substantially the same in our two epochs. 
However, in our sample, 9 of 21 spot-like light curves, observed 
in both epochs, changed their morphology from 2008 to 2011 
and became aperiodic, as seen in Table 5. At the same time, 2 
aperiodic light curves in 2008 became spot-like in 2011. In a 
study of the Orion nebula cluster (ONC), Rebull (2001) analyzed 
Cousins Ic light curves from young stars and showed that in a 
one-year interval, about half of the observed periodic spot-like 
light curves became aperiodic. They also compared their mea¬ 
sured periods with the results of similar analyses of the ONC 
by other authors (Stassun et al. 1999; Herbst et al. 2000) and 
showed that only about 50% of the periodic systems are typically 
recovered between the different works. The results above show 
that cold spots in young stars can evolve substantially and that 
the observable characteristic features, such as light-curve peri¬ 
odicity, can disappear on a timescale of a few years. Cool spots 
are expected to be present in all young low-mass stars, but only 
some of the CTTSs present periodic spot-like light curves, since 
overlying accretion hot spots and circumstellar dust obscuration 
can effectively mask cool spot signatures. When CTTSs are ob¬ 
served to switch between variability categories, it is generally 
attributable to changes in accretion behavior on timescales of a 
few years or less. In Fig. 5 (c) and (d), we present the 2008 and 
2011 light curves of the CTTS Mon-000765, whose light-curve 
classification changed from aperiodic to spot-like. 

Figure 5 shows that flux variations can also be observed from 
2008 to 2011 in the spot-like stars that presented a variable light- 
curve periodicity. During the periodic epoch, most spot-like sys¬ 
tems were fainter than in the aperiodic cycles, as seen in Mon- 
000765 (Fig. 5, lower panels). The brighter phase during the ir¬ 
regular cycles could be due to accretion bursts, leading to the 
appearance of randomly distributed bright spots at the stellar sur¬ 
face, which would at the same time explain the irregular photo¬ 
metric variability. The AA Tau-like stars also show flux variabil¬ 
ity between the periodic and aperiodic phases, but in the oppo¬ 
site sense as compared to the spot-like systems: they are gener¬ 
ally brighter when periodic, as seen in Mon-000928 (Fig. 5, top 
panels). These systems are observed at high inclinations with re¬ 
spect to our line of sight (i > 60°, see McGinnis et al. 2015), and 
the aperiodic cycles may correspond to the unstable accretion 
regime described by Kurosawa et al. (2008) and Kurosawa & Ro¬ 
manova (2013), where multiple accretion streams are randomly 
formed and occult the star as the system rotates. In this scenario, 
the star can be partially occulted at all rotational phases, never 
exhibiting its maximum flux to the observer, and being therefore 
fainter than in the stable and ordered AA Tau phase. This is con¬ 
sistent with the fact that, in the aperiodic cycles, we generally do 
not observe a well-defined maximum in the light curves, while 
in the periodic AA Tau phase a clear maximum is seen, which 
should correspond to the unocculted photospheric brightness. 
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Table 5. Variations of the morphological classification of the 84 CTTS 
light curves observed with CoRoT in 2008 and 2011. 


Light curve 2008 ^ Light curve 2011 N® of stars 


Spot ^ AA Tau 0 

Spot ^ non-periodic 9 

AA Tau ^ spot 0 

AA Tau ^ non-periodic 5 

Non-periodic ^ A A Tau 8 

Non-periodic ^ spot 2 

Non-periodic ^ periodic unclassified 1 

AA Tau ^ AA Tau 8 

Spot ^ spot 12 

Non-periodic ^ non-periodic 37 

Binary ^ binary 2 


Mon-000928 2008 Mon-000928 2011 

15.0 
15.1 
S’ 15.2 

o 

S 15.4 

15.5 

15.6 

4530453545404545455045554560 5890 5900 5910 5920 5930 5940 

JD-2450000 JD-2450000 




(a) 


(b) 


Mon-000765 2008 



JD-2450000 


Mon-000765 2011 

14.04 

14.06 

S’ 14.08 

^ 14.10 

I 14.12 
o 

14.14 
14.16 
14.18 

5890 5900 5910 5920 5930 5940 
JD-2450000 



(C) 


(d) 


Fig. 5. Examples of CTTSs whose morphological classification of the 
light curve varied from 2008 to 2011. Mon-000928 displays an AA- 
Tau-like light curve in 2008 (a) with a period of P = 4.96 days, while 
in 2011 the light curve was aperiodic (b). Mon-000765 was aperiodic 
in 2008 (c), while in 2011 the light curve was spot-like with a period of 
P = 2.41 days (d). 


4.2. Period distribution 

Table 4 shows that ~ 43% of the 96 CTTSs observed in 2008 
and ~ 30% of the 145 CTTSs observed in 2011 were classified 
as periodic. We only classified as periodic the light curves that 
presented a high-power, detached, and dominant peak in the peri- 
odogram analysis. We established a minimum power level of the 
periodogram for a reliable detection based on the periodogram 
analysis of the sinusoidal spot-like light curves, which are clearly 
periodic with a well-determined period within our dataset. We 
also inspected by eye each of the periodic light curves folded 
in phase at the measured period and verified that they presented 
phase coherence. 
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Many light curves that show some periodic signal, but with 
low power in the periodogram, were conservatively classified 
as non-periodic in the present work. Venuti et al. (in prepara¬ 
tion) discuss the peridicity of all the stars observed by CoRoT 
in more detail and the number of periodic stars may be some¬ 
what different from ours. We present in Fig. 6 the period distri¬ 
bution of AA Tau-like and spot-like systems observed in 2008 
and 2011. Despite the observed morphological variations in the 
light curves from 2008 to 2011 (Table 5), the overall period dis¬ 
tribution did not vary significantly from 2008 to 2011. As seen 
in Fig. 6, there are more AA Tau-like stars with long photomet¬ 
ric periods in the new campaign. This is probably due to the 
the difference in datasets between the two campaigns and to the 
longer observation time span in 2011 (40 days) than in 2008 (23 
days), which makes it easier to confidently determine longer pe¬ 
riods in the irregular, though periodic, A A Tau-like light curves. 
Despite the longer observational time span, we did not detect 
longer spot-like periods in 2011. In 2008 and 2011 the spot-like 
light curves were very regular and their periods were easily mea¬ 
sured with our data. Therefore, the lack of periods longer than 
14 days in 2008 and 10 days in 2011 is probably real among 
the observed CTTSs. The longer periods measured were differ¬ 
ent in both campaigns, and this is due to different samples, since 
the stars observed in each run were not exactly the same, and it 
is due to accretion variability, since some spot-like systems be¬ 
came aperiodic and the opposite also occurred between the two 
observational runs. 
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Fig. 6. Period distribution of the CTTSs observed in 2008 (black) and 
2011 (blue). Top: AA Tau-like light curves. Bottom: Spot-like light 
curves. AA Tau-like stars clearly present periods within the range of 
spot-like systems. 


The period of a spot-like light curve is the rotational period 
of the star, since spots are located on the stellar surface. In both 
epochs, the A A Tau-like stars presented periods within the range 
of periods measured in the spot-like light curves. Therefore, the 
material that obscures the AA Tau-like stars is probably not lo¬ 
cated very far from the corotation radius of the star-disk system, 
as pointed out by Alencar et al. (2010). A more detailed anal¬ 
ysis of the period distribution of CTTSs and WTTSs in NGC 
2264 and the implication for angular momentum evolution will 
be presented in another paper (Venuti et al., in preparation). 


5. UV excess and photometric mass-accretion rates 

In this section, we compare the CoRoT light-curve morphol¬ 
ogy with accretion diagnostics and disk indicators to determine 
whether the optical light-curve morphology and accretion and 
inner dusty disk evolution are related. 

Venuti et al. (2014) obtained the UV excess emission of 
CTTSs and WTTSs in NGC 2264 using u- and r -band measure¬ 
ments from CFHT data. They determined the locus of WTTSs 
in the r vs. u-r color-magnitude diagram that was fitted with 
a polynomial function, which was used as a reference value to 
non-accreting systems. The UV excess of a CTTS was defined as 
E(u) = (w-r)obs -(^“Oref^ where (w-r)obs is the observed CTTS 
color, and (u - r)ref is the reference color at magnitude robs- The 
UV excess measurements have an average rms error of 0.16 mag 
to account for the distribution of WTTS around the reference se¬ 
quence (Venuti et al. 2014). With the UV excess, Venuti et al. 
(2014) calculated mass-accretion rates (Muv) for the CTTSs in 
NGC 2264. We show in Fig. 7 the UV excess (top) and mass- 
accretion rate (bottom) vs. Ha equivalent width for the stars in 
our sample. Stars presenting accretion burst light curves show 
some of the highest UV excesses (more negative on the scale 
of the plot) and consequently the highest mass-accretion rates 
among the observed systems, as noted by Stauffer et al. (2014). 
This is supported by a two-sided Kolmogorov-Smirnov test that 
shows that the mass-accretion rate distribution of accretion burst 
stars, when compared to spot-like and A A Tau-like systems, has 
only a 5.8 X 10“^ probability to be equivalent. AA Tau-like and 
spot-like stars are an intermediate population between the accre¬ 
tion bursts systems and WTTSs, which present low UV excess. 
The mean UV excess values of each light curve group corrobo¬ 
rate with this analysis: -1.15 for the accretion bursts, -0.04 for 
the AA Tau systems, 0.05 for the spot-like systems, and 0.14 for 
the WTTSs. These results suggest that the CoRoT light curve 
morphology is related to the evolution of the accretion process. 

Many of the stars observed by CoRoT have Spitzer observa¬ 
tions from the literature (Teixeira et al. 2012), which allowed us 
to compute the aiRAc index, as proposed by Lada et al. (2006). 
This index is the slope of the spectral energy distribution be¬ 
tween 3.6yum and 8/im, and it allows a classification of in¬ 
ner disk evolution. If aiRAc < -2.56, the disk is classified as 
a naked photosphere (no dust in the inner disk). For systems 
with -2.56 < aiRAc < -1.80, the inner disk is optically thin 
and called an anemic disk. For stars with -1.80 < ajRAc < 
-0.5, the inner disk is classified as optically thick. Stars with 
-0.5 < aiRAc < 0.5 are classified as fiat spectrum sources, and 
CHIRAC > 0.5 characterizes class I sources. 

In Fig. 8, we show the aiRAc index vs. the UV excess calcu¬ 
lated by Venuti et al. (2014). Most of the CTTSs present dust in 
the inner disk. Spot-like systems are present in all the inner disk 
classes, but are mostly found among the anemic disks and naked 
photosphere populations. All A A Tau-like systems have dust in 
the inner disk, which is important, since we assumed their light- 
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Fig. 7. Comparison between the CoRoT light-curve classification and 
UV excess (top) or mass-accretion rates (bottom) caleulated by Venuti 
et al. (2014). In the top plot, more negative values indicate higher UV 
excess. 


curve variability to be dominated by circumstellar dust extinc¬ 
tion from that region. Most of the accretion burst systems present 
thick disks, which is consistent with an earlier disk evolution and 
accretion phase than the A A Tau-like stars. The mean ajRAc val¬ 
ues of each light-curve group corroborate this analysis: -1.42 for 
the accretion bursts, -1.71 for the AA Tau systems, -1.84 for 
the spot-like systems, and -2.76, for the WTTSs. These results 
show that the CoRoT light-curve morphology is related to the 
evolution of the inner disk. It also shows that accretion and disk 
evolution are related to each other, since the stars with highest 
UV excesses (more negative values) are the ones with the dusti¬ 
est inner disks, while stars with little or no sign of accretion have 
almost no dust in their inner disks. 

Figure 8 also shows some interesting objects that deserve fur¬ 
ther investigation, which beyond the scope of the present paper, 
however, such as the WTTSs among the anemic and thick disk 
systems. These are non-accreting stars (based on spectra) that 
still show substantial dust in their inner disks. In spite of this, for 
some as yet unknown reason, accretion has been shut off or is 
at a very low level in these systems, at least at the time of spec¬ 
troscopic observation. There is also a significant number of ane¬ 
mic disk candidates among the CTTSs, with various light-curve 
morphologies that would be interesting targets to investigate the 
inner disk dispersal in accreting systems. 
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Fig. 8. Comparison between the CoRoT light-curve classification with 
the UV excess from Venuti et al. (2014) and the slope of the spectral 
energy distribution from 3.6 pm to 8 pm (aiRAc)- These results show that 
the CoRoT light-curve morphology is related to the accretion process 
and the evolution of dust in the inner disk. 

6. Equivalent width and accretion rates from the Her 
line 

The Her line of a CTTS is generally quite strong and wide and 
is often characterized by its equivalent width. Compared to the 
morphology of the CoRoT light curves, the accretion burst sys¬ 
tems show the highest EWh^ in our sample, with an average 
value of 96 + 19 A, clearly distinguishable from the AA Tau and 
spot-like systems. AA Tau-like stars present an average EWh^ of 
13.2 + 2.4 A, while CTTSs with spot-like light curves present a 
mean EWh^ value of 21.1 + 5.6 A. The uncertainties correspond 
to the standard deviations of the measured equivalent width val¬ 
ues. The distribution of EWh^ among the three CTTS groups is 
shown in Fig. 9, where we included EWh^ measurements from 
Dahm & Simon (2005) for the stars that we did not observe. The 
same tendency is observed with the Htr width at 10% of maxi¬ 
mum line intensity. 

Mass-accretion rates (Mhqt) are frequently calculated from 
the Hof line flux, using calibrations such as the ones proposed 
by Fang et al. (2009) and Gullbring et al. (1998). We compared 
the mass-accretion rates obtained from the UV excess by Venuti 
et al. (2014) with the values obtained with the Htr flux, as de¬ 
scribed below. 

We calculated the Htr flux (Fhot) from the EWh^, 


Class I 



lot Spectrum 


Thick Disk ^ 
Anemic Disk 


-O Cold Spot , 
A AA Tau “ 

□ Accretion Bursts 
The other CTTSs 
WTTS 




Naked Photosphere 


Fh. = EWH„Fe(Ha), (1) 

where Fc(Hq') is the continuum flux in the Htr region. Since 
the FLAMES spectra were not flux calibrated, we estimated 
Fc(Hq') from the CoRoT flux using CoRoT observations ob¬ 
tained simultaneously with the FLAMES spectra. The CoRoT 
flux was converted into R magnitudes using the R-band photom¬ 
etry zero-point of 26.74, as determined by Cody et al. (2014). 
For the stars observed with FLAMES that were not observed 
with CoRoT in 2011, we used the continuum flux from the 
CFHT r band (Venuti et al. 2014) taken simultaneously with the 
FLAMES observations. The continuum flux was then corrected 
for extinction with Ay obtained from the CFHT data (Venuti 
et al. 2014). 

The Ha luminosity was then computed as Lh^ = 47rd^FHa, 
where we used d = 760 pc (Sung et al. 1997) as the distance 
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Fig. 10. Comparison between the CoRoT light-curve classification and 
mass-accretion rates obtained from Her flux as a function of Her equiva¬ 
lent width (FLAMESA^LT data). The error bars represent the night-to¬ 
night variability of the line flux. 


Fig. 9. Distribution of Her equivalent width of CTTSs that present AA 
Tau-like (top), spot-like (middle), and accretion burst (bottom) light 
curves. 


to the cluster. We used the fit proposed by Fang et al. (2009) to 
calculate the accretion luminosity (Lace) f^om Lh^: 


log (^) = (2.27 + 0.23) + (1.25 + 0.07) log (2) 

With the accretion luminosity, we obtained the mass- 
accretion rate, 


Mhqt = 


L 


acc 


R* 


GM, 



( 3 ) 


where R* and are the radius and mass of the star, respec¬ 
tively, taken from Venuti et al. (2014), and G denotes the gravi¬ 
tational constant. The inner radius was set to Rin = 5 R=,, which 
is a standard value of the disk truncation radius of CTTSs (Gull- 
bring et al. 1998). The individual values of mass-accretion rates 
calculated from the Her line flux are presented in Table 6. 

Figure 10 shows the mass-accretion rates obtained from the 
mean Her equivalent width of the four to six FLAMES and 
CoRoT simultaneous observations, and also those obtained from 
simultaneous CFHT observations. We plot the mean EWh^, and 
the error bar represents the night-to-night variability of EWh^ 
and Mhot- As observed with the Myv relation (Fig. 7), stars 
that present accretion burst light curves have higher Mh^. Since 
WTTSs are not actually accreting, their Her equivalent widths 
are probably of chromospheric origin, a contribution that is also 
present in the CTTS spectra. Therefore, the mass-accretion rate 
locus of WTTSs represents a lower limit to our ability to measure 
mass-accretion rates of CTTSs from the Her equivalent width in 
our sample, as discussed by Ingleby et al. (2011) and Manara 
etal. (2013). 

Figure 11 shows the relation of the mass-accretion rates cal¬ 
culated with the UV excess and Her equivalent width. The two 
methods show similar tendencies, but the individual value can 
sometimes be different. This difference may arise because the 


Her line does not necessarily originate in the accretion columns 
alone, but can also have contributions from the disk wind and ac¬ 
cretion shock, as absorptions superposed on the accretion profile. 
When we calculated the Her equivalent width, the absorptions 
decrease the measured EWh^ value, leading to lower values of 
mass-accretion rates. At the same time, the errors in Ay affect 
the UV excess accretion rates more than the Her ones and can 
lead to more uncertain values of mass-accretion rates calculated 
from UV excess. Although UV excess is clearly a more direct 
measurement of the mass-accretion rate, Her line fluxes are much 
easier to obtain and allow for a good estimate of mass-accretion 
rates for ah accreting objects, even those that do not present a 
substantial UV excess, due, for example, to circumstehar dust 
obscuration of the hot spot. 



Fig. 11. Comparison between mass-accretion rates obtained from Her 
flux with the UV excess mass-accretion rates. The red line shows a slope 
equal to 1. 


7. Her line variability 

The CTTSs in NGC 2264 present a variety of Her line profiles. 
In this section, we verify how common each Her profile type is 
and compare its occurrence with the CoRoT light-curve classi¬ 
fication. We followed the Her profile classification proposed by 
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Reipurth et al. (1996), which is divided into seven profile types. 
A type I profile is symmetric and shows no absorption feature. 
A type II profile presents two peaks, and the less intense one ex¬ 
ceeds half the intensity of the main peak. A type III profile also 
shows two peaks, and the less intense one is smaller than half the 
intensity of the main peak. A type IV profile shows a pronounced 
absorption feature on the blue or red side of the line where no 
emission is seen, like a P Cygni or inverse P Cygni profile. The 
letters B or R are attached to types II, III and IV, depending on 
the location (in the blue (B) or red (R) side of the line) of the sec¬ 
ondary peak or absorption feature with respect to the main peak. 
We have 20 to 22 nights of observation of each target, and the 
Ha line is quite variable both in intensity and shape from night 
to night. Therefore, we morphologically classified only the mean 
Ha profile of each star and present this classification in Table 6 
and Fig. 12. The top panel of Fig. 12 shows the profile distribu¬ 
tion of all CTTSs observed with FLAMES. Type I is the most 
common Ha profile in our sample. This profile corresponds to 
a main emission feature without any prominent absorption and 
is attributed to magneto spheric accretion. The emission profiles 
of type I stars sometimes display some superposed absorption 
in one or a few observing nights, but the analyzed mean profile 
effectively does not present absorptions. 

Figure 12 (bottom) shows the distribution of the Ha profile 
types according to the CoRoT light-curve classification. Spot¬ 
like stars typically present type I Ha profiles, generally showing 
no evidence of blueshifted absorption due to wind, or redshifted 
absorption due to hot spot photons absorbed by the accretion 
funnel. Moreover, A A Tau-like stars tend to present type HR and 
IIIR profiles, which are characterized by redshifted absorptions, 
compared to the most common type I profile seen among the 
CTTSs in our sample. 

Kurosawa et al. (2006) used disk-wind and magnetosphere 
radiative transfer models for CTTS to examine the Ha classifi¬ 
cation proposed by Reipurth et al. (1996) and with a range of 
parameters reproduced all the Ha profile types. In general, they 
found that the type I profile is the most common among CTTSs, 
since it is produced from a range of system inclinations and is 
dominated by magnetospheric emission. Kurosawa et al. (2006) 
found that type IIB and HR profiles were equally common, as 
also found by Reipurth et al. (1996) and by us, since they re¬ 
sult from the same physical conditions, with the exception of the 
system inclination. Type IIB needs medium inclinations to be 
produced, while type HR results from high inclinations. This is 
also consistent with our results, since most AA Tau-like stars, 
which are high-inclination systems (McGinnis et al. 2015), tend 
to present type HR profiles (see Fig. 12, bottom). 

Reipurth et al. (1996) instead found that the Ha profile dis¬ 
tribution in their sample was dominated by type IIIB profiles 
(33 % of their sample), while in our sample only ~ 9 % of 
the observed CTTSs present type IIIB profiles. In the models 
of Kurosawa et al. (2006), only high mass-accretion rate stars 
presented type IIIB profiles. The difference in profile distribu¬ 
tion found by us and Reipurth et al. (1996) might therefore 
simply be due to a sample selection effect. The stars classified 
as type IIIB by Reipurth et al. (1996) have spectral types be¬ 
tween GO and K7 and a very high Ha equivalent width (average 
value of 64 A), which may indeed indicate that they are high 
mass-accretion rate systems. The type IIIB CTTSs in our sam¬ 
ple present an average EW of 47 A and a mean mass-accretion 
rate of ~ 2.5 x 10“^ Mo/yr, which is not particularly high, but 
does correspond to one of the highest mean mass-accretion rates 
of the various profile types. 
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Fig. 12. Distribution of Her profile types according to the classification 
of Reipurth et al. (1996). Top: All the CTTSs observed with FLAMES. 
Bottom: CTTSs observed both with FLAMES and CoRoT. Light curves 
classified as spot-like, AA Tau-like, and accretion burst are shown as a 
function of profile type. 


8. Periodograms of Her emission iine 

As the star-disk system rotates, the observed Ha line of CTTSs 
varies in intensity and shape. This is expected to occur, even in a 
stable accretion regime, if the stellar magnetic and rotation axis 
are not aligned. The misalignment between the axes breaks the 
axisymmetry of the accretion environment and creates individual 
accretion funnels in each hemisphere. The projection of the Ha 
emission region in our line of sight will then be different at every 
rotational phase, which will cause line variability (see for exam¬ 
ple Alencar et al. 2012). If most of the Ha line is formed in the 
accretion columns, we would expect to measure a periodic line 
variability when the star is in a stable accretion regime. However, 
the Ha line may also form in different circumstellar regions, such 
as the chromosphere, disk wind, and jet, and each region may 
have a different periodicity. This would make a single period 
difficult to detect. With 20 to 22 nights of observation of each 
target with the FLAMES spectrograph, we can check whether 
the Ha line variations were periodic, and compare these periods 
with the periodicities measured photometrically. This will allow 
us to infer if the Ha line-forming region is directly related to the 
physical process causing the photometric variability. 

We analyzed the line periodicity by calculating Scargle pe¬ 
riodograms, as modified by Home & Baliunas (1986), along the 
line profile. The Scargle periodograms fit the data with sinu- 
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soidal functions through a least-squares method. The code re¬ 
quires a range of periods to be tested, and the best fit corresponds 
to the highest power in the periodogram. To obtain the period¬ 
icity in the Ha line, we divided each normalized Ha line profile 
into small velocity intervals (0.5 km/s). For each velocity inter¬ 
val, Scargle periodograms were independently calculated. Each 
periodogram was grouped side by side in velocity, to form a dia¬ 
gram of period vs. velocity. The power of the periodogram is rep¬ 
resented by colors, varying from zero (white) to the maximum 
power intensity (black). With this procedure, we can check how 
each period spreads across the Ha line. Periodic signals only in 
a few velocity bins may not be significant, but periodic signals 
that are spread across several velocity channels indicate period¬ 
icity in an extended region of line formation. We can also verify 
if there are different periodicities in different line regions (like in 
emission and absorption). 

Only eight CTTSs observed with the FLAMES spectrograph 
presented a periodic signal in the Ha line (see Table 6), and their 
periodograms are shown in Figs. 13, 14, and 15. These figures 
are discussed individually in the following subsections, where 
we show that most of the periodicity in the Ha line occurs in the 
red wing, which comes from gas moving away from the observer 
and probably is related to the accretion process. These periodic 
stars presented organized and stable accretion regions rotation- 
ally modulated by the stellar magnetic field. 

Our spectroscopic target selection in 2011 was largely based 
on the CoRoT light-curve classification from 2008, and we in¬ 
cluded as many A A Tau-like and spot-like systems as possible 
in our two pointings to study periodic systems. We did not an¬ 
ticipate, however, that a large number of systems would change 
their light-curve variability from periodic to aperiodic between 
the two runs. Of the 58 stars observed with FLAMES, only 16 
presented periodic CoRoT light curves in 2011, 8 of which were 
classified as A A Tau-like and the other 8 as spot-like systems. 
Probably because of our unfortunate selection, the Ha line of 
most of the 58 observed CTTSs showed no periodicity in our 
data. 

In the next subsections, we analyze the Ha periodicity com¬ 
pared to the CoRoT light-curve morphology discussed in Sect. 4. 
Because of the large number of stars observed spectroscopically 
with FLAMES that presented a non-periodic CoRoT light curve 
in 2011, our main conclusion is that the stars that presented non¬ 
periodic CoRoT light curves are generally not periodic in Ha 
either. We show in Sect. 10 that the lack of periodicity in Ha 
and in the photometry may be caused by an unstable accretion 
regime, as proposed by Kurosawa & Romanova (2013). 

8.1. Ha periodicity in the AA Tau-like stars 

AA Tau-like stars are expected to be in a stable accretion regime, 
where a major funnel flow is formed in each hemisphere, and the 
base is associated with an inner disk warp that periodically oc¬ 
cults the star as the system rotates (Bouvier et al. 2007a; Kuro¬ 
sawa et al. 2008). The photometric period of AA Tau-like stars 
corresponds to the Keplerian period of the occulting circumstel- 
lar material, and it can be different from the stellar rotation pe¬ 
riod, unless the structure responsible for the occultation of the 
star is at the disk corotation radius. As shown in Sect. 4.2, at 
least in our sample, the photometric periods of AA Tau-like sys¬ 
tems is not expected to differ substantially from the stellar rota¬ 
tion period, indicating that the inner disk warp is indeed located 
close to the disk corotation radius. Eight A A Tau-like stars were 
observed with the FLAMES spectrograph, but only three are pe¬ 
riodic in Ha (see Fig. 13). This is somewhat unexpected. If the 


eight stars were in a stable accretion regime, it would naturally 
produce periodic variability in Ha if the line were mostly formed 
in the accretion funnel. However, if the accretion flow is signif¬ 
icantly time variable, it could mask the periodicity imposed by 
the stellar rotation. We will explore evidence for the time vari¬ 
ability of accretion during funnel-flow accretion in a forthcom¬ 
ing paper (Stauffer et al., in preparation). 

The star Mon-000660 (V354 Mon), Fig. 13abc, is an AA 
Tau-like analog that presents the same period in Ha and the pho¬ 
tometry. Assuming that Ha is formed mainly in the accretion 
funnel and that the measured period corresponds to the stellar 
rotational period, it confirms that the occulting disk structure, 
which causes the photometric variability, is located at the coro¬ 
tation radius, as in A A Tau itself (Bouvier et al. 2007a). Another 
AA Tau-like star, Mon-000250, Fig. 13def, is also periodic in Ha 
and presents the same period in the photometry. In the Ha peri¬ 
odogram, this star shows two periodic signals in the red wing, 
the stronger power peak coincident with the photometric period 
(~ 8.9 days) and the weaker peak around ~ 7.5 days. This sec¬ 
ond period is also detected in the CoRoT light curve, but with 
a lower power than the main periodic signal, and is due to the 
periodicity of a small structure present inside the deepest dips of 
the light curve (see Fig. 13f). 

On the other hand, star Mon-000811, Fig. 13ghi, which is 
also AA Tau-like, only shows a faint hint of periodicity in Ha at 
the photometric period (~ 7.9 days), and presents a strong peri¬ 
odic detection at about 13 days in the Ha blue wing. This star 
also presents two periodic signals in the Hel 6678 A line, one at 
the photometric period and another at about 11 days, as shown 
in Sect. 9. The Hel line is expected to be formed near the hot 
spot, close to the stellar photosphere, so one of its periodicities is 
probably close to the stellar rotation period. If the approximately 
eight-day period of Hel corresponds to the stellar rotational pe¬ 
riod, the structure that occults the star could be located close to 
the corotation radius of the star-disk system. This is the only star 
that presents a periodicity only in the blue wing of the Ha line. 
Mon-000811 can be compared to the CTTS SU Aur, which ex¬ 
hibited a blueshifted absorption in Ha that is periodic (Johns & 
Basri 1994), and thought to be due to a wind (Giampapa et al. 
1993). Its variability is explained as rotational modulation of the 
wind by the magnetic field of SU Aur. Mon-000811 shows no 
absorption in the blue wing of Ha (as seen in Fig. 13h), and we 
cannot easily associate the observed periodicity in the blue wing 
with a wind. However, the blue-wing period is longer than the 
Hel period (which we associate with the rotation period), so the 
material responsible for the Ha blue-wing periodicity is proba¬ 
bly located outside the corotation radius and could indeed come 
from a disk wind. 

As discussed above, three AA Tau-like stars present the same 
periodicity in the photometric and spectroscopic analysis. Since 
the spectroscopic periods most likely correspond to the stellar ro¬ 
tational periods, we can assume that, at least in these three cases, 
the inner disk warp is located at the corotation radius. Following 
Bessolaz et al. (2008), we computed the star-disk corotation ra¬ 
dius and calculated the magnetic field strength needed to enforce 
disk truncation at the corotation radius. The results are presented 
in Table 7. At a few stellar radii from the star, the star-disk inter¬ 
action is expected to occur mostly through the dipole component 
of the stellar magnetic field, and the calculated magnetic fields in 
Table 7 are within the range of values of the few dipole compo¬ 
nents measured in CTTSs (B^ip = 0.02 kG to 1.9 kG) (see Table 
2 from Gregory et al. 2012). 
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Table 6. Periods of the classical T Tauri stars observed with the FLAMES spectrograph obtained with photometric and spectroscopic diagnostics. 


Mon ID« 

LC_2008'’ 

PcoRotJ20QS‘’ 

(days) 

LC_2011'’ 

PcoRoT-20n‘’ 

(days) 

Pua" 

(days) 

Puel^ 

(days) 

Regime^ 

Hu Type^ 

(Mo /yr) 

CSIMon-001099 

1 

3.31 

1 

3.38 

NP 

- 

S 

IIB 

3.23e-09 

CSIMon-000103 

1 

1.67 

1 

3.35 

NP 

- 


I 

9.61e-09 

CSIMon-000177 

1 

3.01 

1 

3.01 

NP 

- 


I 

9.29e-09 

CSIMon-000335 

1 

4.51 

1 

4.58 

NP 

- 


I 

2.31e-09 

CSIMon-000965 

1 

9.38 

1 

9.68 

NP 

9.4 

S 

HR 

1.76e-09 

CSIMon-000810 

1 

2.93 

1 

2.93 

NP 

- 


I 

1.24e-09 

CSIMon-000964 

1 

3.34 

1 

3.32 

NP 

- 


I 

1.49e-09 

CSIMon-000326 

1 

7.05 

3 


NP 

NP 


I 

1.99e-09 

CSIMon-000220 

1 

0.76 

3 


NP 

- 


I 

2.50e-10 

CSIMon-000804 

1 

3.23 

3 


NP 

- 


I 

4.83e-10 

CSIMon-001033 

1 

14.15 

- 


6.30 

- 

S 

HR 

4.69e-09 

CSIMon-000250 

2 

4.16 

2 

8.93 

7.5/8.9 

8.6 

S 

HR 

1.44e-09 

CSIMon-000660 

2 

5.25 

2 

5.10 

5.2 

- 

S 

IIIR 

1.45e-10 

CSIMon-000498 

2 

4.23 

2 

4.28 

NP 

- 


I 

l.OOe-08 

CSIMon-001199 

2 

3.75 

2 

3.61 

NP 

- 


IIB 

1.46e-09 

CSIMon-000297 

2 

3.16 

3(E) 


NP 

- 


HR 

5.75e-09 

CSIMon-000928 

2 

4.96 

3 


NP 

10.3 


I 

6.91e-10 

CSIMon-000654 

2 

4.66 

3(E) 


NP 

NP 


I 

4.82e-09 

CSIMon-000441 

2 

4.06 

3(E) 


NP 

Noise 


IIB 

3.24e-09 

CSIMon-000824 

2 

7.05 

- 


NP 

- 

U 

IVR 

3.28e-09 

CSIMon-001144 

3(E) 


3(E) 


NP 

Noise 


II 

2.17e-09 

CSIMon-001275 

3 


3 


NP 

- 


IIIB 

2.51e-08 

CSIMon-000119 

3 


3 


NP 

NP 


TUB 

3.34e-09 

CSIMon-000558 

3(A) 


3 

11.70? 

10.5 

10.7 


IIB 

6.06e-08 

CSIMon-000168 

3(E) 


3(E) 


NP 

Noise 


II 

3.77e-09 

CSIMon-000926 

3 


3 


NP 

NP 


I 

8.66e-09 

CSIMon-000681 

3(E) 


3(E) 


NP 

NP 


IIB 

4.49e-08 

CSIMon-000280 

3 


3 


NP 

- 


HR 

5.90e-10 

CSIMon-000290 

3 


3 


NP 

8.3/11.4 


I 

2.25e-08 

CSIMon-000328 

3 


3 


NP 

NP 


I 

3.90e-09 

CSIMon-000448 

3 


3 


NP 

Noise 


I 

8.16e-09 

CSIMon-000314 

3(E) 


3(E) 


NP 

Noise 


I 

1.06e-08 

CSIMon-000951 

3 


3 


NP 

Noise 


I 

3.21e-08 

CSIMon-000937 

3 


3 


NP 

- 


I 

7.24e-10 

CSIMon-000667 

3(E) 


3(E) 


5.9 

- 

S 

IIB 

1.84e-09 

CSIMon-000945 

3(A) 


3(A) 


NP 

7.4 


I 

6.35e-09 

CSIMon-001022 

3(A) 


3(A) 


NP 

8.9/12.5 


TUB 

6.31e-09 

CSIMon-000996 

3(A) 


3(A) 


8.3 

8.3 


I 

2.79e-08 

CSIMon-000510 

3(A) 


3(A) 


NP 

NP 


I 

2.87e-08 

CSIMon-000341 

3(A) 


3(A) 


NP 

NP 


I 

3.55e-08 

CSIMon-000370 

3 


3 

11.82? 

NP 

NP 


I 

2.55e-08 

CSIMon-000765 

3 


1 

2.41 

NP 

- 


I 

6.48e-09 

CSIMon-000379 

3(E) 


2 

3.68 

NP 

- 

S 

HR 

1.26e-09 

CSIMon-001054 

3(E) 


2 

4.08/8.17* 

NP 

8.1 


I 

1.35e-08 

CSIMon-000296 

3(E) 


2 

3.91 

NP 

- 

S 

HR 

1.73e-09 

CSIMon-000811 

3(E) 


2 

7.88 

12.5 

10.5 

S 

HR 

6.75e-09 

CSIMon-000846 

3 


- 


NP 

NP 


I 

- 

CSIMon-000795 

3 


- 


NP 

NP 


IIB 

6.59e-08 

CSIMon-000457 

3 


- 


NP 

- 


IIIB 

6.01e-09 

CSIMon-000260 

3(E) 


- 


NP 

9.2 


IIB 

2.26e-08 

CSIMon-000893 

3 


- 


NP 

NP 


I 

1.02e-08 

CSIMon-000131 

- 


3 


NP 

9.1 


HR 

- 

CSIMon-000423 

- 


3 


NP 

Noise 


I 

1.44e-08 

CSIMon-000632 

- 


- 


NP 

Noise 


IIB 

1.33e-08 

CSIMon-001287 

- 


- 


NP 

Noise 


HR 

9.38e-08 

CSIMon-000603 

- 


- 


NP 

Noise 


TUB 

6.07e-08 

CSIMon-000239 

- 


- 


5.6 

- 


7 

- 

CSIMon-001128 

- 


- 


NP 

Noise 


I 

1.88e-08 


Notes. This table is ordered by the classification of the CoRoT light-curve morphology from the 2008 campaign. 

CSIMon is an internal identification of the CSI2264 campaign. CoRoT light-curve morphology and photometric period obtained in this work: 1= spot-like, 2= 
A A Tau-like, 3= non-periodic light curves. Accretion bursts stars and aperiodic extinction light curves are identified by (A) and (E), respectively. Period obtained 
with the Hu line from FLAMES spectra. NP denotes a non-periodic line. Period obtained with Hel line 6678A from FLAMES spectra. Noise means that the 
Hel line is present on the spectrum but it is dominated by noise. Accretion regime: S-stable and U-unstable, as discussed in Sect. 10. Morphology of Hu line 
by Reipurth et al. (1996) discussed in Sect. 7. Mass-accretion rates calculated with Hu equivalent width. The photometric period was obtained using Scargle 
periodogram and auto-correlation function, as explained in Sect. 9.1. 


8.2. Ua periodicity in stars classified as non-periodic with the lar variability, but not associated with a main physical process) 
CoRoT photometry showed no well-defined periodicity in their CoRoT light curves 

The stars Mon-000667 (random occultation by circumstellar ma¬ 
terial), Mon-000996 (accretion burst), and Mon-000558 (irregu- 
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Fig. 13. AA Tau-like stars in 2011 that present a periodicity in the Her line. Left column: Bidimensional periodograms of the Her line. The 
color code represents the power of periodogram, varying from zero (white) to the maximum power intensity (black). The black horizontal lines 
correspond to the photometric period of the 2011 CoRoT light curves. Middle column: Her line profiles. Different colors correspond to different 
observation nights, and the thick red line is the average line profile. Right column: CoRoT light curves from 2011. 


Table 7. Parameters used to calculate corotation radius (Rc) and magnetic field strength (B*) of the three AA Tau-like stars that present observa¬ 
tional evidence that the inner disk warp is located at the corotation radius. 


Mon ID« 

M.* 

(Mo) 

RA 

(Ro) 

Mh,." 

(Mo /yr) 

Pcoi?or_201L 

(days) 

(days) 

Rc^ 

(RA 

(kG) 

CSIMon-000250 

1.63 

1.35 

1.44e-09 

8.93 

7.5/8.9 

12.3 

0.81 

CSIMon-000660 

1.86 

1.4 

1.45e-10 

5.10 

5.2 

7.5 

0.07 

CSIMon-000811 

0.91 

1.97 

6.75e-09 

7.88 

7.9/12.5* 

8.2 

0.44 


Notes. 

CSIMon is an internal identification of the CSI 2264 campaign. Stellar mass obtained by Venuti et al. (2014). Stellar radius obtained by Venuti et al. (2014). 
Mass-accretion rates calculated with Hu equivalent width. Photometric period obtained from the 2011 CoRoT light curve. ^ Period obtained with the Hu 
line from FLAMES spectra. Calculated corotation radius. Calculated stellar magnetic field. This star shows only a faint hint of periodicity in Hu at the 
photometric period (~ 7.9 days) that we used as the stellar rotational period. 


in 2008 and in 2011, as shown in Table 6. These stars are peri- Mon-000667 has a period of 5.92 days measured by Lamm 
odic in Ha, however, as seen in Figs. 14a, 14d, and 14g. et al. (2004), which is very similar to our detected period (~ 

5.9 days) and could indeed correspond to the stellar rotational 
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period. In this case, the accreting gas is apparently organized, 
but the inner dust environment does not show a clear correlation 
with the gas structure. 

We could not find any information about the rotational pe¬ 
riod of Mon-000996 in the literature, but it presents an Ha pe¬ 
riod (~ 8.3 days) that is equal to the period found in the Hel 
6678 A line (see also Fig. 18o; we discuss this in Sect. 9). This 
is quite unexpected, since its light curve is dominated by random 
hot spot variability, which would suggest unstable accretion and 
not a globally organized accretion environment. This star also 
presents a similar periodicity (~ 8.1 days) in the Spitzer (IRAC) 
light curve at 3.6yum and 4.5yum (Cody et al. 2014). Recently, 
Blinova et al. (2015) divided the unstable accretion regime into 
two subclasses: ordered unstable (one or two unstable tongues) 
and chaotic unstable (several unstable tongues). A star in the or¬ 
dered unstable accretion regime may retain some periodicity in 
its accretion bursts, as may be the case for Mon-000996, which 
presents only a low-power period detection in its light-curve pe- 
riodogram at about 8.7 days. 

Mon-000558 was classified by Cody et al. (2014), using the 
same CoRoT data, as quasi-periodic, since its light curve in¬ 
cludes low-amplitude stochastic variability superposed on some 
slightly periodic signal. The possible photometric period is ~ 
11.7 days. In the literature, it also presents a period of 0.88 days 
(Flaccomio et al. 2006), which is very different from the Ha and 
the photometric period from our survey, however. We show in 
the next subsection that the Ha period of ~ 10.5 days is also 
seen in the Hel line (see Table 6), which may indicate that it cor¬ 
responds to the stellar rotation period. This star had an accretion 
burst CoRoT light curve in 2008, while in 2011 it has an irreg¬ 
ular, possibly hot-spot light curve. This is another system that 
is a candidate to be in the ordered unstable regime described by 
Blinova et al. (2015). 

8.3. Ha periodicity in stars not observed by CoRoT in 2011 

Two stars that show Ha line periodicity were not observed by 
CoRoT in 2011. The star Mon-001033, Fig. 15a, although not 
observed by CoRoT in 2011, shows a possible period of 14 days 
measured in 2011 with the Spitzer (IRAC) light curve in the 
CSI2264 campaign (Cody et al. 2014). This star was observed 
by CoRoT in 2008, and its light curve was classified as spot¬ 
like, with a period of 14.15 days, which agrees with the IRAC 
period, but is more than twice that obtained with the Ha line 
(~ 6.30 days). The periodic region in Ha corresponds to a nar¬ 
row velocity range, where a redshifted absorption is present in 
the spectrum. 

The other star that presented a periodicity in Ha, but was 
not observed by CoRoT, is Mon-000239. This star does not have 
any period information in the literature and is not periodic with 
Spitzer photometry. The detected period in Ha occurs in a very 
limited profile region associated with a small redshifted absorp¬ 
tion feature, like the one observed in Mon-001033. 

8.4. Ha periodicity in spot-like stars 

Eight CTTSs classified as spot-like with the 2011 CoRoT light 
curves were observed with FLAMES. However, despite being 
periodic in the photometry, none of them showed any periodicity 
in the Ha line. Two examples are presented in Fig. 16. The pho¬ 
tometric periodicity of spot-like stars is mostly due to cold spots 
on the photosphere, while the Ha line is mainly formed in the 
accretion funnel. Results from spectropolarimetry have shown 


that the main hot and cold spots tend to coincide in accreting 
stars (Donati et al. 2010, 2011a,b), and we would therefore ex¬ 
pect the cold spot and the main accretion funnel, which ends in 
the hot spot, to present the same periodicity. At the same time, 
the Ha line profiles of spot-like systems tend to be weak and 
may therefore have a significant chromospheric contribution. In 
some spot-like systems, the lack of Ha line periodicity could 
then be due to a rather uniform chromospheric emission con¬ 
tributing significantly to the profile, which would dilute the Ha 
modulation from the accretion funnel. 

9. Periodograms of the HeI6678A emission iine 

The Hel 6678 A singlet line is often found in emission in stars 
with a moderate to high accretion rate (Muzerolle et al. 1998; 
Beristain et al. 2001b; Luhman et al. 2003). Like the Hel 5876 A 
triplet line, it may present a narrow (NC) and a broad (BC) com¬ 
ponent, with different kinematic characteristics. In the spectra of 
CTTSs, the NC is found to be slightly redshifted and circularly 
polarized and is interpreted to arise from gas in the post-shock 
region near the stellar surface (Beristain et al. 2001b; Edwards 
et al. 2003). The BC is more complex. It can be found to be red¬ 
shifted, but is most commonly observed blueshifted in CTTSs. 
The current interpretation is that it may come from the accre¬ 
tion funnel or the base of a hot accretion-powered stellar wind 
(Beristain et al. 2001b; Edwards et al. 2006). 

Depending on the star-disk inclination with respect to our 
line of sight, the hot spot may be hidden by the inner disk and 
remain invisible most of the time. In that case, we would not 
expect to see the Hel line in emission, even at high accretion 
rates. This is probably the case of the A A Tau-like system Mon- 
000660, which is seen at high inclination (79° + 11°, Fonseca 
et al. 2015) and does not present UV excess or Hel 6678 A in 
emission, despite being actively accreting, as shown by its ex¬ 
tended Ha emission profile. 

The Hel NC is formed close to the stellar surface, therefore 
its periodicity, whenever present, should correspond closely to 
the stellar rotation period. We expect to detect periodic variabil¬ 
ity if accretion is in a stable regime, which creates a main accre¬ 
tion funnel and consequently a main hot spot in each hemisphere. 

We found that 24 of the 58 CTTSs observed with FLAMES 
clearly presented the Hel 6678 A line in emission. Another set of 
11 CTTSs also show Hel in emission, but the S/N is too low in 
the line to allow a reliable analysis of line profile variability. Ac¬ 
cording to Beristain et al. (2001b), the Hel NC is characterized 
by line widths smaller than 50km/s. In our sample, the width 
of the mean Hel 6678 A line profile ranged from 29 to 47 km/s, 
depending on the target, with a mean value of 39.3 km/s. There¬ 
fore, in our sample, all of the stars that showed Hel 6678 A 
in emission only presented an NC. Close to the Hel 6678 A 
line, there is a photospheric absorption Eel line (at 6677 A). We 
used residual spectra to calculate the Hel 6678 A line equivalent 
widths because of this line. The Eel 6677 A line was subtracted 
from these spectra, using as template the spectrum of a WTTS 
of the same spectral type of each star. 

We calculated the Hel line periodograms as was done for 
the Ha line (see Sect. 8). The periodograms of the Hel lines are 
shown in Figs. 17 and 18. Of the 24 CTTSs with clear Hel emis¬ 
sion, 12 stars showed no periodicity in this line, and they were 
also non-periodic in the CoRoT photometry. The non-periodicity 
of half of the Hel emission stars may indicate the intermittent 
formation of several short-lived hot spots near the stellar surface, 
instead of one well-organized hot spot in each hemisphere. 


Article number, page 14 of 33page.33 



A. R Sousa et aL: CSI2264: Accretion process in classical T Tauri stars in the young cluster NGC 2264 


Mon-000667 



(a) 


Mon-000996 



(d) 


Mon-000667 



(b) 


Mon-000996 



v(km/s) 

(e) 


Mon-000667 2011 



JD-2450000 

(C) 


Mon-000996 2011 



JD-2450000 

(f) 




-400 -200 0 200 400 

V (km/s) 


(g) 


Mon-000558 



(h) 


Mon-000558 2011 


13.25 

CJt 

o 
2 
CH 

I 13.30 
o 
o 

13.35 

5890 5900 5910 5920 5930 5940 
JD-2450000 

(i) 



Fig. 14. Periodic stars in the Her line that are not periodic in the CoRoT photometry. Left column: Bidimensional periodograms of the Her line. 
The color code represents the power of periodogram, varying from zero (white) to the maximum power intensity (black). Middle column: Her line 
profiles. Different colors correspond to different nights of observation, and the thick red line is the average line profile. Right column: CoRoT light 
curves from the 2011 campaign. 


Of the 12 systems that showed a periodicity in the Hel line, 
two presented more than one period in this line (as seen in Table 
6 and Figs. 18e and 18i), which may be because of the formation 
of several long-lived (compared to the time span of our observa¬ 
tions) hot spots at different latitudes on the stellar surface, rather 
than a main hot spot on each hemisphere. In the next subsections, 
we present the analysis of the Hel line periodicity compared to 
the CoRoT light-curve morphology discussed in Sect. 4. 

9.1. Hel line in AA Tau like stars 

AA Tau-like stars are seen at high inclination with respect to 
our line of sight, and the inner disk is expected to hide part 
of the stellar photosphere as the system rotates. Depending on 
the system geometry, the hot spot, and consequently the region 
where the Hel NC forms, may also be hidden in these systems. 
Only three of eight A A Tau-like stars observed with FLAMES 


showed the Hel line in emission, and these three systems pre¬ 
sented a periodicity in this line. The Hel line periodograms of 
Mon-000250, Mon-001054, and Mon-000811 are shown in Fig. 
17 and the measured periods are presented in Table 6. 

Mon-000250 shows the same periodicity in Hel and Ha lines 
(Figs. 17a and 13d), and also presents the same period in the 
CoRoT light curve, which confirms that the structure responsi¬ 
ble for the stellar occultation is located close to the co-rotation 
radius of the star-disk system. 

Mon-001054 (Fig. 17c) shows a Hel period of ~ 8.1 days, 
that is, twice the CoRoT photometric period from 2011. The 
photometric period, calculated with the modified Scargle pe¬ 
riodogram, as explained in Sect. 4, presents a low detection 
power. Using a different period search method, such as the auto¬ 
correlation function (Cody et al. 2014; McQuillan et al. 2013), 
we find a ~ 8.2 day period in the CoRoT light curve, but the de¬ 
tection significance is again very low. From the literature, Mon- 
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Fig. 15. Periodic stars in the Her line that were not observed by CoRoT in 2011. Left column: Bidimensional periodograms of the Her line. The 
color code represents the power of the periodogram, varying from zero (white) to the maximum power intensity (black). Middle column: Her line 
profiles. Different colors correspond to different nights of observation, and the thick red line is the average line profile. Right column: CoRoT light 
curve from 2008, when available. 
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Fig. 16. Spectroscopic and photometric variability of spot-like systems. Left column: Bidimensional periodogram of the Her line of two spot-like 
systems. The color code represents the power of periodogram, varying from zero (white) to the maximum power intensity (black). Black lines 
correspond to the photometric period of the stellar light curve observed with CoRoT in 2011. Middle column: Her line profiles. Different colors 
correspond to different nights of observation, and the thick red line is the average line profile. Right column: CoRoT light curve from 2011. None 
of the spot-like systems observed with FLAMES showed any periodicity in the Her line. 
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Fig. 17. Periodic stars in the Hel 6678 A that also showed a periodicity in the CoRoT light curves. Bidimensional periodograms of the Hel line 
(first and third columns). The color code represents the power of the periodogram, varying from zero (white) to the maximum power intensity 
(black). The black horizontal lines correspond to the period of the 2011 CoRoT light curve. The second and fourth columns present the Hel 6678 A 
lines used to calculate the periodograms. Different colors correspond to different nights of observation, and the thick red line is the average line 
profile. 


001054 has a photometric period of ~ 7.8 days (Lamm et al. 
2004), which agrees with our Hel period and could correspond 
to the stellar rotation period. 

Mon-000811 (Fig. 17e) is periodic in Hel (~ 10.5, higher 
power and ~ 8 days lower power) and Her (~ 12.5 days) (Figs. 
17e, 13g) and shows periodicity in the CoRoT light curve (~ 
7.88 days). However, these periods are different from each other. 
As discussed above, this star is classified as AA Tau-like, and 
we would therefore expect accretion to occur in a stable regime, 
if the inner disk warp is associated with the base of the accre¬ 
tion column, and forming a major hot spot in each hemisphere 
at the top of the accretion column. One possible scenario to ex¬ 
plain this complex system is that the Hel line eight-day period 
corresponds to the stellar rotation period, which is close to the 
photometric period, and therefore the inner disk warp that causes 
the photometric variability is located close to the co-rotation ra¬ 
dius. The Hof ~ 12.5-day period is observed in the blue wing of 
the emission profiles and could then correspond to emission or 
absorption from a wind located outside the corotation radius. 

9.2. Hel line periodicity in spot-like stars 

Because of the high stability of the spot-like light curves on the 
timescale of our observations and the low mass-accretion rates 
of these systems, we initially assumed that large cold spots were 
the main cause of the observed photometric variability of spot- 
like systems. If, instead, hot spots were causing the photometric 
variability, we might expect to observe the same type of vari¬ 
ability in the Hel line. We observed eight spot-like systems with 


FLAMES, and only one, Mon-000965, has the Hel line in emis¬ 
sion, as seen in Table 6. The Hel periodogram of this star is 
shown in Fig. 17g, and it presents a clear periodicity that co¬ 
incides with the period obtained with the CoRoT photometry. 
Mon-000965 is not clearly periodic in Ha (Sect. 8.4), but a faint 
periodic signal (low power in the periodogram) is seen at the 
photometric period. The fact that this star has the Hel NC line in 
emission indicates that the hot spot is prominent, which suggests 
that its light-curve variability might be dominated by a stable hot 
spot instead of a cold spot, as initially assumed during our 2011 
campaign. 

9.3. Periodic stars in the Hel line that are not periodic in the 
CoRoT photometry or were not observed by CoRoT in 
2011 

Half of the 24 stars that showed Hel in emission presented pe¬ 
riodicity in this line, but only 4 (Sect. 9.1 and 9.2) were also 
periodic in the CoRoT photometry, as seen in Fig. 17. The other 
8 stars that showed a periodicity in the Hel line were not found to 
be periodic in the CoRoT light curves of 2011 (6 stars) or were 
not observed with CoRoT in 2011 (2 stars), as seen in Fig. 18. 

Ah of the accretion burst systems observed by FLAMES 
(five in total) are included in the stars that showed Hel in emis¬ 
sion, which is reasonable, since they ah presented high mass- 
accretion rates and had light curves dominated by hot spot vari¬ 
ability (Stauffer et al. 2014). Three of these five stars (Mon- 
000996, Mon-001022 and Mon-000945) presented a periodicity 
in the Hel line, despite showing no periodic signal in the CoRoT 
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Fig. 18. Periodic stars in the Hel 6678 A line that did not show a periodicity in the CoRoT light curves. The first and third columns show 
bidimensional periodograms of the Hel line. The color code represents the power of the periodogram, varying from zero (white) to the maximum 
power intensity (black). The second and fourth columns show the Hel 6678 A lines used to calculate the periodograms. Different colors correspond 
to different nights of observation, and the thick red line is the average line profile. 


photometry. Mon-000996 (18o) also presented similar periodic- 10. Stable and unstable accretion 
ities in both the Hel and Ha lines, as shown in 14d. 
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CTTSs are photometrically variable on different timescales, 
from seconds to even decades (Rucinski et al. 2008), depend¬ 
ing on the various physical phenomena that affect the stellar and 
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circumstellar environment. Variable accretion is one of the main 
sources of variability in this evolutionary phase. As discussed 
by Kurosawa & Romanova (2013), accretion can occur in sta¬ 
ble and unstable regimes. The unstable regime may be triggered 
by the Rayleigh-Taylor instability, which takes place in the inner 
disk region. This type of variability produces non-periodic light 
curves that are dominated by accretion bursts, as observed in 
some stars in NGC 2264 (see Stauffer et al. 2014). This variabil¬ 
ity also produces detectable spectral features in emission lines, 
such as HyS, Hy, and PayS (Kurosawa & Romanova 2013). 

According to theoretical predictions, stars in a stable accre¬ 
tion regime will present redshifted absorption in just a few rota¬ 
tional phases, when the main hot spot is seen projected through 
the accretion funnel in our line of sight. In this situation, photons 
emitted by the hot spot will be absorbed by the high-velocity gas 
in the accretion funnel, causing the observed redshifted absorp¬ 
tion components. In the unstable accretion regime, redshifted ab¬ 
sorption may be seen at any rotational phase because of the large 
number of accretion funnels, and consequently, hot spots on the 
stellar surface. The accretion burst stars are all expected to be 
in the unstable accretion regime. Unfortunately, only five accre¬ 
tion burst systems were included in our sample of stars observed 
with FLAMES, and none of them showed redshifted absorption 
in Ha. Their Ha line profiles are dominated by emission due to 
accretion (Reipurth type I profile) that typically does not present 
redshifted absorption (see Fig. 12). To investigate the occurrence 
of redshifted absorption in the emission lines of accretion burst 
systems, we need to analyze higher order Baimer lines, as sug¬ 
gested by Kurosawa & Romanova (2013), or the Hel (10830 A) 
line, which often shows both redshifted and blueshifted absorp¬ 
tion in CTTS spectra, as discussed by Edwards et al. (2006) and 
Cauley & Johns-Krull (2014). 

In the simulations by Kurosawa & Romanova (2013), the 
redshifted absorption components are not visible in the Ha line 
either because its emission is very intense and hides the shallow 
absorption. Redshifted absorption is most commonly seen ob- 
servationally in weaker Baimer emission lines such as H^ and 
Hy. We do not have the higher B aimer lines in our spectra, but 
ten stars show a clear redshifted absorption component in the Ha 
line (marked in Col. 8 of Table 6). Despite the small number of 
targets, we tried to verify how the redshifted absorption compo¬ 
nent varies with stellar rotational phase and compare our results 
with the prediction of stable and unstable accretion simulations. 

Of our ten selected targets, only Mon-000824 is possibly in 
the unstable accretion regime, where the redshifted absorption 
is present in the line profile at all rotational phases. This star 
was classified as AA Tau-like in 2008, but was not observed by 
CoRoT in 2011 (see Fig. 19), which means that in principle we 
would not know its photometric behavior when the spectroscopic 
data were obtained. As seen in Sect. 4, CTTSs are highly variable 
and in a few years may change the morphology of the their light 
curve from AA Tau to non-periodic and vice versa. Although 
Mon-000824 was not observed with CoRoT in 2011, during the 
CSI campaign about 900 epochs of I-band photometry were ob¬ 
tained from November 2011 to March 2012 of the central region 
of NGC 2264, including Mon-000824, with the US Naval 40” 
telescope. As discussed by McGinnis et al. (2015), the simulta¬ 
neous I-band and CoRoT photometry typically match very well. 
We therefore inspected the I-band light curve of Mon-000824 in 
2011 and found that it was indeed not periodic, implying that 
this star probably changed from a stable to an unstable accretion 
regime from 2008 to 2011. 

The other nine selected stars (two spot-like, five AA Tau- 
like, one aperiodic extinction variables and one not observed 



-200 0 200 -200 0 200 -200 0 200 -200 0 200 
v(kin/s) 


Fig. 19. Ha line of Mon-000824 as a function of rotational phase, in¬ 
dicated in each plot. This star was not observed by CoRoT in 2011, 
and we used the photometric period (7.05 days) obtained in 2008, when 
the star presented an AA Tau-like light curve, to calculate the rotational 
phases. This star is a candidate to be in the unstable accretion regime in 
2011 . 


by CoRoT in 2011) are apparently in a stable accretion regime, 
showing redshifted absorptions in only a few rotational phases. 
In Fig. 20, we show the spectra of Mon-000296 as an example. 

The predictions that stable and unstable accretion regimes 
present a different rotational phase distribution of redshifted ab¬ 
sorption components are based on the theoretical line profiles 
computed by Kurosawa & Romanova (2013) for a system with 
an inclination of 60° with respect to the line of sight. However, in 
systems with lower inclinations, hot spots will always be in view 
(Kurosawa et al. 2008), hence, even in a stable accretion regime 
the redshifted absorption may be present in the line profile at all 
rotational phases. This is, however, not the case of Mon-000824 
(Fig. 19), since its inclination with respect to our line of sight is 
72° + 5° (McGinnis et al. 2015). 

We would expect stars in a stable accretion regime to show 
periodic light curves and stars in an unstable accretion regime to 
show non-periodic light curves (Kurosawa & Romanova 2013; 
Kulkami & Romanova 2008, 2009). This is because the stable 
accretion regime generates an organized circumstellar environ¬ 
ment that shows periodicity as the star-disk system rotates. Of 
the nine stars we classified as being in a stable regime, using this 
criterion of redshifted absorption, only two did not present pe¬ 
riodic light curves, as seen in Table 6. Instead their light curves 
were classified as aperiodic extinction variables. The stable ac¬ 
cretion regime classification therefore agrees with the observa¬ 
tion of periodicity in the CoRoT light curves for most stars with 
redshifted absorption in just a few rotational phases. 

In the stable accretion regime, the variability of the red¬ 
shifted absorption is expected to be periodic. Of the nine stars 
we classified as being in a stable regime, five are periodic in the 
Ha line, agreeing with the theoretical predictions of Kurosawa 
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Fig. 20. Hof line of Mon-000296 as a function of rotational phase, in¬ 
dicated in each plot. The central peak is the nebular line that was not 
entirely removed. This star presents an AA Tau-like light curve with a 
3.91 day period and is a candidate to be in the stable accretion regime 
in 2011. 

& Romanova (2013). This means that only about half of the stars 
we classified as being in a stable accretion regime are periodic in 
the Ha line. The only star in the unstable accretion regime (Fig. 
19) is also not periodic in the Ha line, as predicted by Kurosawa 
& Romanova (2013). 

The classification of a stable accretion regime based on 
the presence of redshifted absorption at well-defined rotational 
phases agrees only sometimes with the type of observed vari¬ 
ability in the CoRoT light curves. Neither is a star that shows 
redshifted absorption at all rotational phases necessarily accret¬ 
ing in the unstable regime. Other factors, such as the inclination 
of the system, may influence the presence of the redshifted ab¬ 
sorption in the Baimer lines. Moreover, we must keep in mind 
that the Ha line is not the best choice for this type of analysis, 
since other lines, such as H/3, Hy, and FajS, are weaker, and the 
redshifted absorption is more commonly visible in the line pro¬ 
files. 

11. Correlation matrices 

The Ha emission line can be formed in several different regions, 
in the accretion funnel, the chromosphere, and/or the wind. The 
analysis of the correlation of different parts of the Ha line profile 
can help us investigate if the line comes from different regions 
and if these regions are correlated. 

One method of checking the correlation across a line profile 
is to calculate correlation matrices. We divided each Ha profile 
into small velocity intervals of 1.5 km/s and calculated the linear 
correlation coefficient, r(ij), of the line intensity at each velocity 
interval i with each velocity interval j, as described by Johns & 
Basri (1995a). The correlation coefficient ranges from -1 to 1, 


1 corresponding to a perfect correlation, 0 to no correlation, and 
-1 to a perfect anticorrelation. When i = y, riij) = 1, and for 
all values of i and y, r(/,y) = r(y,/), implying that the matrix 
is symmetrical relative to the diagonal. We present the result as 
a color-map plot of the correlation coefficients, which is useful 
to visually identify regions of the profile that are correlated or 
anticorrelated. This has been used in the literature to analyze 
correlations across the Ha line and other lines such as Hp (Johns 
& Basri 1995b; Oliveira et al. 2000; Alencar & Batalha 2002; 
Kurosawa et al. 2005). 

All of the 58 CTTSs that have FLAMES spectra present 
some correlation signal of the red wing with itself. Some stars 
show this correlation just in the redshifted absorption position 
and others just in the redshifted emission. We find that ~ 74% 
of the CTTSs show a very good correlation across the entire Ha 
red wing, which then varies coherently. Almost all the CTTSs 
observed by FLAMES also exhibited some correlation of the 
blue wing with itself, ~ 74% presenting a substantial correla¬ 
tion. The blue wing of Ha also varies coherently. This shows 
that the variability in each wing of the Ha profile is coherent, 
which indicates that each line wing (red or blue) is formed in a 
region dominated by a given physical process. 

We find that ~ 34% of the CTTSs observed with FLAMES 
show no sign of correlation between the red and the blue wings, 
implying that they vary independently in those systems. At the 
same time, ~ 66% CTTSs show some correlation of the red wing 
with the blue wing and among these, -55% show substantial 
correlation across the entire wings, while - 45% show correla¬ 
tions that occur only in some specific velocity range. 

Only 7 of the 58 CTTSs observed with FLAMES present an¬ 
ticorrelated variations in the Ha line profile. The anticorrelation 
seen in 4 of these 6 stars is associated with a blue- or redshifted 
absorption that varies in antiphase with the rest of the profile. 
The other three stars present an anticorrelation associated with a 
blueshifted emission and the other profile parts. 

Overall, the Ha line profiles present variations that cannot 
be explained by continuum variability alone that is due, for ex¬ 
ample, to hot spots. If excess continuum emission from hot spots 
were the main cause of emission line variability, the entire profile 
would be correlated, and the Ha line would not vary in shape, but 
only in intensity. Other factors such as wind emission or absorp¬ 
tion, accretion variability, and non-symmetrically distributed cir- 
cumstellar material must influence the emission line variability, 
often simultaneously, since many stars do not present a very 
good correlation across the entire Ha line profile. It is possible, 
however, that in some cases the Ha line may be formed mainly 
in one specific region, like the wind or the accretion column. We 
would then expect the line profile to be dominated by the same 
physical process and the variability of the entire line profile to 
be correlated. 

Figures 21 and 22 show examples of correlation matrices 
for some CTTSs, together with their observed Ha line profiles. 
These figures are representative of the different types of corre¬ 
lation matrices we obtained. In these plots, the correlation coef¬ 
ficients are represented by colors ranging from black (r = -1, 
maximum anticorrelation) to light blue (r = 0, no correlation) 
and to orange (r = 1, maximum correlation). The correlation 
matrices and Ha line profiles of all the CTTSs in our sample 
observed with FLAMES are presented in the Appendix. 

Figure 21a,b shows that for Mon-000457 the red wing of Ha 
correlates well with itself in the velocity range of 50 km/s < v < 
350 km/s, and the blue wing shows a good correlation with itself 
at -350 km/s < v < -75 km/s. We can also see a region of cor¬ 
related variability in the blue wing at -100 km/s < v < Okm/s, 
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which corresponds to the blueshifted absorption from a wind 
that varies quite independently of the rest of the profile. The 
blue wing is anticorrelated with the red wing for velocities in 

the range of-100 to-350km/s and ~ 0 to ~ 350km/s. 

The line profiles show that these regions correspond to the blue 
and red emissions that vary in antiphase with each other. This 
indicates a strong influence of different physical processes in 
each Ha profile wing, the red wing is probably dominated by 
accretion and the blue wing by the wind. The accretion process 
and the wind may contribute to the entire profile, but can have a 
stronger contribution in one particular velocity range, as seen for 
example in theoretical line profiles calculated by Kurosawa et al. 
(2006) and Lima et al. (2010), which include both the wind and 
accretion contributions. In the case of Mon-000457, if the blue 
wing is mostly influenced by the wind and accretion is proba¬ 
bly the main process acting on the red wing, the anticorrelation 
between the two wings could be explained as being due to rota¬ 
tional modulation of a CTTS with an inclined magnetosphere, as 
also observed in SU Aur (Johns & Basri 1994). When the main 
accretion funnel faces the observer, accretion is strong and ejec¬ 
tion is low, and half a rotation cycle later, the opposite occurs. 
In an ideal case, these variations should be periodic, but the Ha 
line variability of this star does not show very clear signs of pe¬ 
riodicity. 

Figure 21c shows a strong correlation of the red wing with 
itself and with the blue wing of the Ha profile at high veloci¬ 
ties (|v| > 200km/s) for Mon-001022 . At the same time, the 
strong blueshifted absorption (-200km/s < v < Okm/s) is an¬ 
ticorrelated with the blue and red emission regions of Ha. This 
feature is due to photons from the accretion columns and a hot 
spot absorbed by the outflowing disk wind material in our line 
of sight. Stars that present very strong blueshifted absorption in 
Ha normally present emission profiles dominated by the wind, 
as can be seen in the theoretical profiles calculated by Kurosawa 
et al. (2006) and Lima et al. (2010). When the emission due to 
the wind increases, photons are even more strongly absorbed by 
the wind, increasing the depth of the blueshifted absorption. This 
causes the blueshifted absorption to become anticorrelated with 
the rest of the emission line. 

Figure 22a shows the correlation matrix and the Ha line pro¬ 
files of Mon-001033. The blue and red wings of the Ha profile 
correlate well with themselves, but show no correlation or an¬ 
ticorrelation with each other. The emission component is quite 
variable, not very intense, and suffers strong influence from ab¬ 
sorption in the blue and red wings. Blueshifted absorption is as¬ 
sociated with winds and redshifted absorption with the accretion 
funnel. The combination of faint emission and strong absorption 
that probably comes from different physical processes created a 
variability pattern where the variations of one wing are not cor¬ 
related with the other wing. 

Theoretical correlation matrices are rare in the literature. 
In Alencar et al. (2012), Ha correlation matrices were calcu¬ 
lated from theoretical line profiles that included only the mag- 
netospheric contribution. Their theoretical correlation matrices 
showed a good correlation of the entire Ha line profile, but out¬ 
flows were not taken into account. Our observed CTTSs include 
stars such as Mon-000945, whose Ha emission line profiles 
show no hint of blueshifted absorption and should be dominated 
by accretion. Figure 22c,d shows the Ha profiles of Mon-000945 
and its corresponding correlation matrix, which indeed is well 
correlated across the entire profile, as predicted theoretically by 
magnetospheric accretion models. However, the Ha profile of 
most CTTSs cannot be explained by magnetospheric accretion 
alone, and it would be very interesting to be able to compare 



Fig. 21. Examples of correlation matrices (left) of the Her line of two 
CTTSs and the corresponding line profiles (right). In panels (a) and (c), 
the color range corresponds to the value of the linear correlation coeffi¬ 
cient between the different velocity bins of the Her line profiles. Perfect 
anticorrelation corresponds to -1 (black), no correlation to 0 (light blue), 
and a perfect correlation to 1 (orange). In panels (b) and (d), different 
colors correspond to spectra observed in different nights, and the thick 
red line is the average line profile. The correlation matrices and Her line 
profiles of all the CTTSs in our sample observed with FLAMES are 
presented in the Appendix. 


the observed correlation matrices with theoretical matrices cal¬ 
culated from profiles that include both accretion and outflows. 

Veiling variations would also produce a correlation across 
the entire line profile, as observed in the Her line of Mon-000945. 
However, the Her line profile presents changes in its morphology, 
not only in its intensity (see Fig. 22d), as would be expected 
for variations that are only due to the extra veiling continuum. 
The type of correlated variability and Her profile morphology ob¬ 
served in Mon-000945 is also seen in the CTTSs T Tau (Johns & 
Basri 1995b) and TW Hya (Alencar & Batalha 2002), for which 
it has also been suggested that the Her line is mostly formed in 
the accretion columns. 

The analysis of the correlation matrices shows that the Her 
line variability can be associated with different physical pro¬ 
cesses. The variability of the red wing often does not affect the 
variability of the blue wing. This is confirmed by the lack of cor¬ 
relation of the blue wing with the red wing, seen in ~ 34% of the 
Her correlation matrices. This lack of correlation disagrees with 
theoretical models that only include the magnetospheric accre¬ 
tion contribution to the emission line profile. It would be very 
interesting to compare observational correlation matrices to ma¬ 
trices calculated from hybrid magnetosphere and wind models, 
such as those of Kurosawa et al. (2006, 2011), and determine 
which circumstellar environment would be able to reproduce the 
observational matrices. At the same time, the red wing is almost 
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Mon-001033 



Mon-000945 



Fig. 22. Examples of correlation matrices (left) of the Her line of two 
CTTSs and the corresponding line profiles (right). The color code is the 
same as in Fig. 21. 


always well correlated with itself, as is the blue wing with itself, 
indicating that the variability of each wing is dominated by the 
same physical phenomenon. 


12. Interesting objects 

Some objects are interesting enough to be analyzed individually 
in future work. We found three spectroscopic binaries with the 
FLAMES data. The CTTS Mon-000448 is a double-line spec¬ 
troscopic binary, while the WTTS Mon-000497 and the CTTS 
Mon-000804 are single-line spectroscopic binaries. These stars 
do not present eclipses in the CoRoT light curves, however. 

The AA Tau-like star Mon-00811 is a CTTS that has differ¬ 
ent periods in different accretion and circumstellar features. One 
possible interpretation is that the different periods come from 
different regions of the star-disk system. It has a photometric pe¬ 
riod (~ 7.9 days) due to occultation of the star by the inner disk, 
a stellar rotation period (~ 8 days) obtained from the Hel line 
periodogram, a 10.5-day period from the Hel blueshifted emis¬ 
sion, and a period of ~ 13 days in the Ha line blue wing that we 
associate with a disk wind. 

Some WTTSs have strong IR excesses, but do not show 
a clear sign of accretion (Mon-001157, Mon-000434, Mon- 
000279, Mon-000753, and Mon-000271). For some unknown 
reason, accretion has been shut off or is at a very low level, de¬ 
spite the available material in their inner disks. It would be in¬ 
teresting to reobserve these systems to verify whether accretion 
features reappear in the future, of if it has indeed come to an end. 
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13. Conclusions 

We investigated accretion in a group of classical T Tauri stars 

that belong to the young cluster NGC 2264, using observational 

data obtained in the Coordinated Synoptic Investigation of NGC 

2264 campaign. The main results of this work are listed below. 

1. The light curves observed by CoRoT were morphologi¬ 
cally classified as spot-like, AA-Tau-like, or non-periodic. 
Of the non-periodic light curves, some are due to accre¬ 
tion bursts, others to circumstellar disk obscuration, but most 
non-periodic light curves are very complex and could not be 
associated with a major physical phenomenon. In the 2011 
campaign, we found ~ 13 % spot-like, ~ 12 % AA Tau-like, 
and ~ 70 % non-periodic systems. Of the non-periodic ones, 
~ 19 % were associated with accretion bursts, ~ 13 % with 
occultation by circumstellar dust, and ~ 68 % presented vari¬ 
ability of unknown origin. 

2. Of the CTTSs observed by CoRoT, ~ 43% and ~ 30% 
showed periodic photometric variability in 2008 and 2011, 
respectively. The number of periodic stars did not change 
significantly between the two runs, but a larger number of 
non-periodic stars were included in the 2011 campaign. 

3. Of the 84 stars that were observed with the CoRoT in the two 
different epochs (2008 and 2011), ~ 30% showed substan¬ 
tial variability in their light-curve morphology, changing to 
a different category in our light-curve classification scheme. 
These changes show how dynamic the CTTS circumstellar 
environment is in only a few years. 

4. When UV excesses are not available, mass-accretion rates 
can be obtained with good results from the Ha flux using the 
currently available calibrations from the literature. 

5. Most of the CTTSs we observed showed no periodicity in 
the Ha line. Only eight stars in our sample of 58 CTTSs 
observed with FLAMES were periodic in Ha. 

6. The Hel 6678 A narrow component line variability can pro¬ 
vide the rotation period of a CTTS when the star is in a stable 
accretion regime, since in that case a major hot spot is ex¬ 
pected to form in each stellar hemisphere. We found that 24 
of the 58 stars observed by FLAMES had this line in emis¬ 
sion. Of these, 12 were periodic in the Hel line, and 4 were 
also periodic photometrically. 

7. None of the photometrically periodic spot-like systems that 
were observed spectroscopically presented a periodicity in 
Ha, while spectropolarimetry results show that hot and cold 
spots tend to be coincident in CTTSs, and we might expect 
to observe the same periodicity in the cold spot and accre¬ 
tion funnel diagnostics. Spot-like systems tend to present low 
mass-accretion rates, and the lack of periodicity in the Ha 
line variability of spot-like systems could be due to chromo¬ 
spheric emission contributing significantly to the emission 
profile in these systems, which would dilute the Ha modula¬ 
tion from the spot. 

8. We analyzed which CTTSs were in stable or unstable ac¬ 
cretion regimes. Kurosawa & Romanova (2013) suggested 
that CTTSs in a stable accretion regime should present red- 
shifted absorption components in emission line profiles only 
in a few rotational phases, when photons from the hot spot 
are seen in our line of sight projected against the main accre¬ 
tion column. Systems in unstable accretion will produce ran¬ 
dom accretion tongues and should present redshifted accre¬ 
tion components at any rotational phase. We analyzed the ten 
CTTSs that showed redshifted absorption components in Ha 
and for which we had some information about the rotational 
period. Of these, nine CTTSs fulfilled the stable accretion 
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criteria (five stars were AA Tau-like, two stars were spot¬ 
like, one was non-periodic, and one was not observed in 2011 
by CoRoT), and in almost all the systems accretion in a sta¬ 
ble regime agrees with the type of variability observed in the 
CoRoT light curves. Only one star of the ten selected CTTSs 
was apparently in the unstable accretion regime, presenting 
redshifted absorption at all rotational phases and showing no 
periodicity in Her, as expected. 

9. The analysis of correlation matrices showed that the variabil¬ 
ity in the Her line can be associated with different physical 
processes, and the red and blue wings are sometimes formed 
in different regions, since the variability of one wing does 
not always correlate with the variability of the other. This 
is confirmed by the lack of correlation of the blue with the 
red wings, seen in ~ 34% of the CTTSs. Each wing tends 
to be well correlated with itself, indicating that each wing is 
dominated by a single physical phenomenon. 
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Table 2. Parameters of the CTTSs observed by CoRoT and/or FLAMES/VLT. 


Mon ID« 

2Mass ID^ 

Spr 

EWna" 

(A) 

ErEw^ 

(A) 

(kms”*) 

EWhc" 

(A) 

FUR06^ 

O^IRAC^ 

(2008) 

LC" 

(days) 

p/i 

(2011) 

LC" 

(days) 

p/i 

Q' 

CSIMon-000007 

06415304-^0958028 

K7 





c 

-1.47 

3(A) 


3(A) 


0.75 

CSIMon-000011 

06411725-^0954323 

K7 




58.3 



3(A) 


3(A) 



CSIMon-000017 

06413199-^1000244 

K5 




13.1 


-2.69 

1 

4.78 

1 

4.73 

-0.03 

CSIMon-000021 

06405944-^0959454 

K5 




7.4 

c 


2 

3.23 

2 

3.16 

0.15 

CSIMon-000056 

06415315-^0950474 

K5 




1.8 

w 

-1.36 

2 

5.71 

2 

5.86 

0.14 

CSIMon-000058 

06420870-^0941212 

K4.5 




94.0 

c 


- 


1 

2.14 

0.22 

CSIMon-000063 

06411193-^0959412 

M2.5 




19.4 


-1.37 

- 


3 



CSIMon-000090 

06410896-^0933460 

M3 




51.0 

c 

-0.88 

3 


3 



CSIMon-000103 

06405954-^0935109 

K6 

15.6 

0.9 

401 

6.4 

c 

-1.14 

1 

1.67 

1 

3.35 

0.21 

CSIMon-000117 

06405413-^0948434 

M2.5 




353.0 


-2.34 

- 


3(A) 


0.72 


Notes. Only a portion of this table is shown here. A full version is available in the online journal. This table is ordered according to the Mon ID. 
(a) “cSIMon” is an internal identification of the CSI 2264 campaign. Elsewhere in the text "CSI" was omitted for brevity. 2MASS identification. Spectral type 
obtained by Venuti et al. (2014); Dahm & Simon (2005); Rebull et al. (2002); Walker (1956). Hu parameters obtained in this work using FLAMES spectra. We 
used the convention that positive Hu equivalent width indicates Hu in emission, and negative values indicate Hu in absorption. The uncertainties were obtained 
assuming a Poisson distribution. Hu equivalent width obtained by Dahm & Simon (2005). ^ Classification as CTTS (c) and WTTS (w) by Furesz et al. (2006). 

Chirac is the slope of the spectral energy distribution between 3.6pm and 8//m obtained by Teixeira et al. (2012). CoRoT light curve morphology and 
photometric period obtained in this work: “1”= spot-like, “2”= AA Tau-like, “3”= non-periodic light curves. Accretion bursts and aperiodic extinction light curves 
are identified by (A) and (E), respectively. “UNC” = unclassified light curve and “B”= Binary. Parameter defined by Cody et al. (2014) to distinguish periodic 
from aperiodic light curves. According to Cody et al. (2014), when Q>0.6 the star is not periodic photometrically. 
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Table 3. Parameters of the WTTSs observed by CoRoT and/or FLAMES/VLT. 


Mon ID« 

2Mass ID" 

SpT^ 

EWHa" 

FrEw^ 

W10%Ha‘' 

EWh„" FUROe' aiRAc^ 




(A) 

(A) 

(kms“') 

(A) 


CSIMon-000008 

06414859+0954115 

K5 



CSIMon-000009 

06420914+0948048 

F5 



CSIMon-000014 

06420664+0941317 

K7:M0 



CSIMon-000015 

06420911+0959027 

K7:M0 



CSIMon-000018 

06411322+0955086 

K3:K4 

3.0 

-2.73 

CSIMon-000020 

06420924+0944034 

K7 


w 

CSIMon-000022 

06410457+1000426 

M4 



CSIMon-000024 

06415684+0947451 

Ml 



CSIMon-000026 

06410518+1000189 

M4 

6.5 


CSIMon-000029 

06410328+0957549 

K7 

1.5 

-2.85 


Notes. Only a portion of this table is shown here. A full version is available in the online journal. This table is ordered aceording to the Mon ID. 
(a) “cSIMon” is an internal identification of the CSI 2264 campaign. 2MASS identification. Spectral type obtained by Venuti et al. (2014); Dahm & Simon 
(2005); Rebull et al. (2002); Walker (1956). Hu parameters obtained in this work using FLAMES spectra. We used the convention that positive Hu equivalent 
width indicates Hu in emission, and negative values indicate Hu in absorption. The uncertainties were calculated assuming a Poisson distribution. Hu equivalent 
width obtained by Dahm & Simon (2005). ^ Classification as CTTS (c) and WTTS (w) by Furesz et al. (2006). (Xirac is the slope of the spectral energy 
distribution between 3.6 jum and 8jum obtained by Teixeira et al. (2012). 
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Appendix A: Correlation matrices and Ha line 
profiles 

We present in this appendix the correlation matrices and Ha line 
profiles of all the CTTSs we observed with the FLAMES spec¬ 
trograph during the CSI 2264 campaign. We organized the cor¬ 
relation matrices according to the three main morphologies dis¬ 
cussed in Sect 11. Some correlations matrices are difficult to fit 
in a single morphological class, and in these cases, we classified 
the star according to the most representative feature present in 
the correlation matrix. In Fig. A.l, we grouped all the CTTSs 
that show anticorrelation in some part of the Ha line profile, like 
the stars in Fig. 21. A few stars that were not included in this 
first group present correlation matrices with anticorrelation only 
in the narrow nebular emission region, near zero velocity. This is 
due to difficulties in removing the nebular contribution, and the 
anticorrelation is not related to an emission or absorption region 
of the star-disk system. In Figs. A.2, A.3 and A.4, we organized 
the stars that show correlation in almost all of the profile, like 
the matrix in Fig. 22c. Finally, in Figs. A.5 and A.6, we present 
the stars that show no sign of correlation between the red and the 
blue wings in most of the profile, like the matrix in Fig. 22a. 
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Fig. A.l. Correlation matrices (left) of the Her line of selected CTTSs and the corresponding line profiles (right). In the left panels, the color range 
corresponds to the value of the linear correlation coefficient, r(/,j), between different velocity bins (/ and j), of the Her line profiles, as described 
in Sect. 11. Perfect anticorrelation corresponds to -1 (black), no correlation to 0 (light blue), and a perfect correlation to 1 (orange). When i = j, 
r(ij) = 1, and for all values of i and j, r(ij) = implying that the matrix is symmetrical relative to the diagonal. In the right panels, 

different colors correspond to spectra observed in different nights, and the thick red line is the average line profile. The stars in this group present 
anticorrelation of some part of the line profile. The anticorrelation seen in these stars is associated with blue- or redshifted absorption or emission 
that varies in antiphase with the rest of the profile. For example, when in panel a (Mon-000457) the red wing emission increases in intensity, the 
blue wing emission decreases. 
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Fig. A.2. Same as Fig. A.l, but for stars that present strongly correlated line profiles, as indicated by the positive correlation in almost the entire 
line profile. The narrow anticorrelation region in the matrices of Mon-000632, Mon-001612, Mon-000119, Mon-000314, and Mon-000448 is due 
to the nebular contribution that was not entirely removed. 
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Fig. A.3. Same as Fig. A.2 
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Fig. A.4. Same as Fig. A.2 
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Fig. A.5. Same as Fig. A.l, but for stars that show no sign of correlation between the red and blue wings of the Ha line profile. The narrow 
anticorrelation in the matrices of Mon-000824, Mon-000177, and Mon-000667 is due to the nebular contribution that was not entirely removed. 
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Fig. A.6. Same as Fig. A.5. The anticorrelation in the matrice of Mon-000250 is due to the nebular contribution that was not entirely removed. 
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Fig. A.7. Same as Fig. A.5. The anticorrelation in the matrice of Mon-000795 is due to the nebular contribution that was not entirely removed. 
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